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ABSTRACT 


This  thesis  develops  and  evaluates  practical  methods  of  'online' 
excitation  control  for  power  system  stabilization. 

The  equations  of  a  power  system  are  expressed  in  a  state  space 
form.  A  method  of  finding  the  closed  loop  quasi -optimal  excitation 
control  is  presented.  The  linearized  model  of  a  single  machine  infinite 
bus  problem  is  considered  and  the  method  is  then  extended  to  a  nonlinear 
model.  The  problem  of  implementing  the  quasi -optimal  control  is  studied 
and  a  practical  compromise  is  suggested.  The  results  obtained  with  the 
quasi -optimal  scheme  are  compared  with  those  obtained  by  a  steepest  descent 
method.  It  is  found  that  the  quasi -optimal  scheme  gives  optimal  solution 
for  linear  systems.  The  effect  of  time  delay  on  the  stability  of  the 
system  is  also  considered.  Finally,  the  closed  loop  (quasi -optimal )  control 
has  been  applied  to  a  two  machine  system  and  the  effect  of  local  control 
investigated. 
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CHAPTER  1 


INTRODUCTION 

The  growth  of  loads  as  well  as  increased  demand  for  reliability  have 
led  to  the  development  of  complex,  highly  interconnected  power  systems. 

» 

Disturbances  on  such  a  system  may  propagate  throughout  the  rest  of  the 
system,  resulting  in  major  failures  of  power  supply  in  absence  of  adequate 
safeguards.  For  reliable  production  and  transmission  of  electrical  energy, 
it  is  essential  that  the  power  system  be  stable  in  the  sense  described  in 
the  following  section.  System  stability  can  be  enhanced  by  suitable 
control  techniques,  some  of  which  are  described  in  Section  (1.2).  One  of 
the  more  important  of  these  —  control  of  machine  field  voltage  or  excit¬ 
ation  control  is  investigated  in  this  thesis.  A  corresponding  strategy 
for  state  feedback  quasi -optimal  excitation  control  is  developed  and 
evaluated  in  the  chapters  to  follow. 

1.1  Power  System  Stability 

The  definition  of  power  system  stability  as  given  by  Kimbark^  is 
"Power  system  stability  is  a  term  applied  to  alternating  current  electric 
power  systems,  denoting  a  condition  in  which  the  various  synchronous 
machines  of  the  system  remain  in  synchronism  or  'in  step'  with  each  other. 
Conversely,  instability  denotes  a  condition  involving  loss  of  synchronism 
or  falling  'out  of  step'".  Depending  on  the  magnitude  and  type  of 
disturbances,  power  system  stability  may  be  classified  into  three  types. 

These  are  described  below. 

Steady  state  stability  is  the  ability  of  the  power  system  to  maintain 
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power  transfer  over  the  system  without  loss  of  stability  or  synchronism 
when  the  magnitude  of  power  transfer  is  increased  gradually.  It  is 
assumed  that  the  increase  in  power  level  upto  this  limit  occurs  slowly 
enough  to  allow  regulating  devices  to  respond  with  their  steady  state 
characteristics,  and  that  inertia  effects  are  negligible. 

Transient  stability  is  the  ability  of  the  power  system  to  maintain 
stability  in  the  presence  of  a  sudden  large  change  in  load  occasioned  by 
system  switching  (as  may  be  caused  by  load  switchings  or  circuit  breaker 
operations),  or  by  a  fault.  It  is  usually  assumed  that  the  regulating 
devices  do  not  have  time  to  function  during  the  transient  period  and  that 
non-linear  modes  of  system  operation  are  encountered  during  the  period. 
Inertia  is  important. 

Dynamic  stability  is  the  ability  of  the  power  system  to  maintain 

stability  for  small  disturbances  and  to  prevent  growth  of  oscillations. 

Dynamic  instability  generally  occurs  due  to  lack  of  damping  torque.  Due 

to  the  small  disturbance  assumption,  this  type  of  stability  is  usually 

[2  3] 

investigated  through  the  use  of  linearized  models1-  * 

1.2  Basic  "Swing  equation11  and  Methods  to  Improve  Stability 

The  "Swing  equation"  for  a  synchronous  machine  can  be  written  as 

"4tDl  =Pin-Po  (1J> 

dt 

where  the  output  power  PQ  is  a  function  of  the  rotor  angular  position, 
machine  internal  voltages  and  transfer  (or  mutual)  admittances  of  the 
system.  The  right  hand  side  of  equation  (1.1)  represents  the  acceler¬ 
ating  power  P  .  At  a  steady  state  operating  point  of  the  machine, 
a 

both  the  velocity  of  the  machine  5  and  the  accelerating  power  Pfl  are  zero 
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while  the  rotor  angle  and  other  machine  currents  are  constant.  This  is 
also  called  an  equilibrium  point.  The  following  approaches  have  been 

[41 

usedL  J  to  help  stabilize  a  system  after  a  perturbation  from  the  equilibrium 
point  takes  place  due  to  a  disturbance. 

1.  Changing  transfer  admittance:  Any  switching  of  the  system 

i 

configuration  to  improve  power  transfer  between  sources  is  beneficial.  The 
faster  this  is  done  the  better,  (a)  The  most  commonly  used  method  involves 
detection  and  clearing  of  faults  by  protective  relays  and  circuit  breakers. 
The  time  required  has  been  decreased  through  the  years  due  to  the  use  of 
improved  equipment,  to  the  order  of  3  cycles  at  the  present  time.  Any 
further  improvement  in  this  figure  is  not  likely,  (b)  The  transfer 
admittance  of  the  power  system  can  be  changed  by  switching  in  additional 
capacitors  into  the  transmission  network.  Control  of  power  system  transients 
by  capacitor  switching  has  been  reported  in  1  iterature^’^  in  recent  years. 
Rama  Rao^  and  Miniesy^  have  considered  the  problem  of  finding  the 
optimal  value  of  capacitance  for  particular  disturbances.  The  difficulty 
encountered  by  these  investigators  has  been  in  pre-determining  the  optimal 
value  of  capacitance  and  the  duration  of  switching  independent  of  the  type 
of  disturbance,  (c)  Braking  resistors  are  also  used  to  change  the 
effective  conductance  of  the  system  and  for  dissipating  the  extra  energy 
available  during  a  fault.  Load  shedding  is  useful  in  decreasing  load 
on  heavily  loaded  generators. 

2.  Change  in  mechanical  input:  Pin  varies  in  accordance  with  the 
speed  governor  action  following  deviations  from  scheduled  frequency  or  as 
a  result  of  supervisory  control  action.  Since  mechanical  systems  have 
large  time  constants,  stabilization  by  controlling  Pin  is  not  very  effective. 


♦  ’ ... 

. 


4 


However,  "fast  valving"  may  be  helpful  in  this  regard^-L 

3.  Change  in  field  excitation:  This  directly  controls  the 
internal  voltage  of  the  machine  thus  controlling  its  power  output.  This 
approach  will  be  considered  in  detail  in  the  next  section  since  this  type 
of  control  for  stabilization  is  used  in  this  thesis. 

1.3  Historical  Development  of  Excitation  Control 

Until  recently  most  synchronous  machines  were  equipped  with 
rotating  type  of  exciters  having  relatively  long  time  constants  and  low 
"ceiling"  voltages.  As  a  result  they  were  of  limited  effectiveness  in 
improving  the  stability  of  the  power  system. 

The  introduction  of  static  excitation  systems  has  revolutionized 
the  technology  of  excitation  control.  In  recent  years,  considerable 
emphasis  has  been  placed  on  the  use  of  static  excitation  schemes  with 
relatively  fast  response  and  high  ceiling  voltages  to  improve  power  system 
stability.  Two  major  effects  of  high  speed  voltage  regulators  are  well 
known y^g  fjrsy  js  the  increase  in  the  restoring  synchronizing 
forces  made  possible  through  the  forcing  of  both  excitation  and  internal - 
machine  fluxes.  This  improves  transient  or  "first  swing"  stability.  The 
second  is  the  deterioration  of  machine  damping  resulting  in  dynamic 
instability.  The  use  of  additional  feedback  signals,  particularly  velocity, 
through  the  excitation  system  has  been  suggested  for  satisfactory 
operation^  .  Most  studies^0,  ^  have  considered  linear  models  even 
when  dealing  with  first  swing  stability.  For  large  disturbances  a  non¬ 
linear  model  is  needed.  It  has  been  found  possible  to  select 
control  signals  such  as  the  magnitude  of  acceleration  that  effectively 
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control  large  disturbances  but  degrade  dynamic  stabi  1  i  ty^ ^ .  Dineley^^ 
found  that  by  a  'continuous  judicious  combination1  of  velocity  and 
acceleration  signals  a  maximum  effect  on  rotor  oscillations  can  be  produced. 
Jones  and  Smith^]  considered  the  problem  of  transient  removal  with 
bang-bang  excitation  control. 

1.4  Research  Objective 

The  equations  governing  the  power  system  dynamics  are  of  high  order 
and  nonlinear.  To  obtain  the  optimal  control  for  stabilization  (the 
performance  index  for  such  problems  is  either  time, or  the  norm  of  the  devia¬ 
tion  of  some  or  all  of  the  states  and  controls  from  the  nominal  value)  for 
such  a  problem  by  standard  optimization  procedures,  one  has  to  resort  to 
iterative  techniques  even  for  disturbances  known  a-priori ^ >8, 16 ,17] , 

However,  disturbances  which  appear  on  a  power  system  are  in  general  not 
known  in  advance.  Some  system  parameters  such  as  the  transfer  admittances 
are  dependent  on  the  location  and  type  of  disturbance,  and  have  to  be 
measured  immediately  after  a  disturbance  appears  on  the  system.  The  control 
must  again  act  within  a  very  small  fraction  of  a  second  to  stop  'first 
swing'  instability  for  large  disturbances.  By  the  time  the  optimal  control 
is  obtained  for  the  particular  system  configuration  by  means  of  'off-line' 
computations,  it  may  be  too  late  and  instability  may  result.  The  time 
limitation  is  not  very  severe  for  small  disturbances,  however,  since  the 
oscillations  grow  slowly. 

Thus,  in  general,  a  truly  optimal  excitation  control  scheme  may 
not  be  implementable  due  to  the  complexity  of  the  controller  and  time 
limitations.  For  this  reason  it  is  considered  worthwhile  to  explore  quasi - 
optimal  excitation  schemes  which  are  implementable  on  real  systems. 
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The  main  objective  of  this  study  is  to  explore  means  for  stability 
enhancement  by  eliminating  system  transients  in  the  least  possible  time. 

For  such  control  schemes  to  be  realizable  'on-line'  it  is  necessary  to 
develop  techniques  for  efficient  determination  of  control  in  terms  of 
system  state  and  parameters.  Secondary  objectives  are  investigation  of 
the  implementability  of  such  control  and  consideration  of  practical 
compromises  needed  to  apply  the  control  to  real  systems. 

1 .5  Scope  of  Thesis 

The  quasi -optimal  excitation  control  for  the  linearized  single 
machine  infinite  bus  problem  is  found  directly  as  a  function  of  the  system 
states.  Pontryagin's  minimum  principle  is  used  to  obtain  the  control 
from  a  transformed  simplified  system.  A  proportional  control  is  suggested 
as  a  practical  compromise  which  for  small  disturbances  is  equivalent  to  a 
velocity  feedback  signal.  The  scheme  is  then  extended  to  a  nonlinear 

/ 

machine  model  and  a  two  machine  system.  Results  indicate  that  the  quasi- 
optima!  excitation  control  is  very  effective  in  removing  transients,  though 
a  bang-bang  excitation  scheme  may  not  be  realizable.  The  proportional 
control  which  is  a  modification  of  the  bang  bang  control  is  realizable  and 
effective.  A  summary  of  results  and  conclusions  as  well  as  suggestions 
for  future  research  will  be  presented  in  Chapter  6  . 
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CHAPTER  2 

MATHEMATICAL  MODEL  OF  A  POWER  SYSTEM 

The  synchronous  generator,  exciter  and  load  models  are  presented 
in  the  first  three  sections  of  this  chapter.  The  mathematical  models  for 
a  power  system  are  derived  in  the  state  space  form  in  Sections  4  and  5. 
Two  cases  are  considered:  (a)  A  single  machine  connected  to  an  infinite 
bus  (hereafter  referred  to  as  single  machine  case),  (b)  two  machines 
feeding  a  load  from  a  busbar. 

2.1  Synchronous  Generator  Equations 


Fig.  2.1  Multi -machine  System  Configuration 


11 


Assume  that  p  machines  are  interconnected  by  a  transmission 

network  so  that  they  feed  a  single  equivalent  load.  A  multi -load  system 

n  21 

can  be  reduced  to  this  form.  The  equations  for  the  mth  machine  areL  *  J 
(The  list  of  symbols  is  given  in  page  xiii) 
a)  Voltage  and  flux  linkage  equations  for  the  armature 
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-i  - 

(2.2) 


*  Subscript  m  refers  to  the  quantities  of  the  mth  machine. 
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b)  Voltage  and  flux  linkage  relations  for  the  rotor 
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Here  by  definition 


xafdo 

rfd 


V 


fd 


(2.5) 


It  is  assumed  that  there  are  j  rotor  windings  on  each  axis  and  a  field 
winding  on  the  direct  axis. 

c)  The  equation  for  the  terminal  voltage  of  the  mth  machine  is 


e 


tm 


(2.6) 


d)  The  system  connection  relationships  for  the  mth  machine  with  the 
bus  bar  are 


(2.7) 


where 


v,  =  v  sin(<$  +  e  ) 
d  m  v 

v  =  v  cos(<$  +  0  ) 
q  m  v 


(2.8) 


is  the  angle  between  bus  bar  voltage  and  a  common  reference  frame. 
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e)  The  torque  equations  for  the  mth  machine  are 


T 

r 

0  1 

• 

^m^m  ^in^rn 

4> 

-1  0 

• 

_q 

m 

_  qj 

(2.9) 

(2.10) 


In  order  to  distinguish  between  reference  frames,  subscripts 
d  and  q  are  used  to  denote  components  along  the  two  axes  of  machine 
reference  frames  while  subscripts  D  and  Q  refer  to  components  of  the 
common  reference  frame.  The  currents  and  voltages  of  the  machines  are 
related  to  the  network  (having  common  reference  frame)  by  the  trans¬ 
formations 


"VD 

COS  6 
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-sinSm 
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-V 
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sin6 

m 

cos  6 

rrr 

V 

COS5 

m 

sin6m 

m 

h. 

-sin6,„ 
u  m 

cosS 

nW 

(2.11) 


(2.12) 


2.2  The  Excitation  System 

There  are  various  types  of  voltage  regulators  in  use  but  only  the 
direct  acting  continuous  type  of  regulator  and  exciter  is  considered  here. 
The  power  amplification  for  modern  continuous  regulators  is  provided  by 
such  devices  as:  magnetic  amplifiers,  amplidynes  and  controlled  rectifiers. 
Since  the  voltage  regulator  controls  the  reactive  power  of  a  generator,  it 


.  . 
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usually  includes  certain  limiting  and  compensating  circuits.  These 
additions  are  not  considered  here.  When  a  linearized  model  of  the  system 
is  used,  the  exciter  and  regulator  may  be  represented  mathematically 
by  a  series  of  gains  and  time  constants^.  In  this  study,  however,  the 
voltage  regulator  and  exciter  is  represented  by  a  single  gain  and  time 
constant  even  for  nonlinear  modes  of  system  operation. 

Static  exciters  or  electronic  exciters  are  now  almost  universally 
used  for  large  generators.  The  maximum  output  voltage  (or  ceiling  voltage) 

f5l 

for  these  exciters  is  as  large  as  7  per  unit1-  J  compared  to  about  2  per 
unit  for  conventional  type  of  exciters.  The  corresponding  time  constants 
are  also  of  the  order  of  10-20  milliseconds  compared  to  about  .5  seconds 
for  that  of  an  amp! i dyne  type  of  exciter.  This  type  of  exciter  may  be  very 
reasonably  represented  by  a  single  gain  and  time  constant.  The  voltage 
regulator  equation  is  written  as 

AEfd  =  TTTJJ  (et  -  etr  '  “.<*»  (2J3) 

where  K  is  less  than  zero, 
r 

2.3  Load  representation 

Various  load  representations  are  used  for  stability  studies. 
Constant  impedance,  constant  current  or  constant  power  and  constant  power 
factor  load  models  may  be  used.  Until  recently  loads  in  stability  studies 
have  been  customarily  represented  as  impedances  which  are  constant  under 
changing  voltage  and  frequency  conditions.  This  practice  has  been  accepted 


r  5  .  s  •'  - 


v  • 
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for  one  or  more  of  the  following  reasons^ 

(a)  For  studies  of  the  first  transient  swing  of  a  small  system 
characterized  by  a  generating  plant  feeding  radially  into  a  load  area,  the 
load  representation  is  not  very  critical. 

(b)  The  constant  impedance  elements  lend  themselves  to  simple 
models.  In  this  case  loads  may  be  included  in  the  admittance  matrix  of 
the  system.  Until  the  advent  of  large  scale  digital  computers,  it  has 
been  impractical  to  make  any  other  assumption  for  load  response. 

(c)  Data  on  the  response  of  system  load  to  disturbances  has 
been  generally  unavailable  in  the  literature  and  it  is  almost  impossible 
to  obtain  such  data  on  the  system. 

(d)  It  is  possible  to  achieve  a  completely  closed  form^ 
simulation  algorithm  by  making  the  assumption  that  all  loads  have  a 
constant  impedance  current-voltage-frequency  characteristic  which  saves 
computational  time  and  effort. 

From  the  above  considerations,  the  simplest  type  of  load 
representation,  constant  impedance,  is  used  throughout  the  study.  The 
load  is  described  by  the  complex  matrix  equation 

[V]  -  [Z][I]  (2.14) 

Or  by  its  inverse 

[I]  =  [Y][V]  (2.15) 

Each  of  these  equations  can  be  expanded  to  a  set  of  the  form 
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(2.16) 


2.4  Single  machine  infinite  bus  system 

In  stability  studies,  the  interval  of  time  of  interest  is  relatively 
small  compared  to  the  time  constants  of  the  governor  and  the  prime  mover. 
While  governor  and  prime  mover  responses  are  important,  they  are 
ineffective  in  the  early  part  of  the  machine  swing  and  will,  if  considered, 
contribute  to  stability.  For  simplicity,  the  torque  variations  due  to 
governor  control  are  neglected.  The  amortisseur  windings  are  known  to  aid 
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stability,  so  these  are  also  not  considered.  The  system  studied  is 
shown  in  figure  2.2. 


Fig.  2.2  System  configuration  for  single  machine  case. 

For  this  case  ev  in  equation  (2.8)  is  zero.  Combining  the  terminal 
constraint  with  the  machine  flux  linkage  equations  and  replacing  by 
u(t),  the  following  set  of  ordinary  nonlinear  differential  equations  is 
obtained.  (Details  are  given  in  Section  A2-1  of  appendix  A2). 


X  =  f(x)  +  d_  u(t) 

where 

*  ■  [ifd’  v  y  n> s] 


d  is  a  vector  of  constant  elements 
The  voltage  regulator  equation  is 


pu(t)  =  y  et  -  f-  -  7"  us(t)  +  (Eo  '  Kr  etr)/xr 
r  r  r 


(2.17) 

(2.18) 


(2.19) 


The  linearized  state  equations  for  the  machine  given  in  section  A2-2  of 
appendix  A2  are 
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AX.  =  £  A_X  +  d_  All  ( t) 


(2.20) 


The  corresponding  linearized  regulator  equation  is 


ACi(t)  =  l 

j=l 


C6j  Xj 


K 

+  C66  Au(t)  -  ^  U  (t) 

r 


(2.21) 


2.5  Two  machine  system 

The  following  system  configuration  is  considered 


machine  n 


Fig.  2.3  System  configuration  for  two  machine  case. 

As  in  the  case  of  the  single  machine,  the  governor  action  and  amortisseur 
windings  are  neglected.  The  following  set  of  differential  equations  is 
obtained 


X  =  f (X  ,  v)  +  d  um(t) 
An  — -m  — m  nr 


(2.22) 
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The  torque  equations,  in  particular,  are 


^in^m  ^xafd  Vd  iq^m  +  ^xd  ”  xq^m  ^d 


q'm 


6m  =  a)  n 
m  o  m 


(2.23) 

(2.24) 


where  m  =  1 ,  2 

Equation  (2.22)  is  similar  to  equation  (A2-1.13)  except  that  a-^  is  a 
function  of  the  variable  bus  voltage  v  and  arguments  of  sine  and  cosine 
terms  change  from  6  to  6+e^.  The  bus  bar  voltage  is  obtained  by  the 
relation 


v 


=  Z  1 


n 


where  i  is  the  load  current. 

& 


(2.25) 


... 
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APPENDIX  A2-1 


The  differential  equations  for  the  single  machine  infinite  bus 
case  are  developed  in  this  section.  For  this  case  ©v  in  equation  (2.8) 
is  zero.  It  is  assumed  that  there  is  only  one  winding  on  the  rotor  -  the 
field  winding. 

Substituting  the  terminal  constraints  equations  (2.7)  and  (2.8) 
in  equation  (2.1)  and  rewriting  equation  (2.3)  gives  (neglect  saturation) 


r  worfd 

p<lfd  '  Efd  '  “o  rfd  nfd 


(A2-1.1) 


pib  ,  =  u)  [—  pi  .  +  R  1  , - XI  +  VSin6] 

K^d  oLo)  rd  e  d  u>  e  q  J 

o  o 


+  Wq  R  id  +  w*d 


(A2-1.2) 


pill  =  01  [—  pi  +  R  i  +  —  x  i  .  +  vcoss] 
oLco  r  q  e  q  to_  e  a 
^  0  0 


+  “0  R  \  ■  “♦d 


(A2-1.3) 


Substituting  the  flux  linkage  relations  (2.2)  and  (2.4)  in  the  above 
three  equations,  the  following  simultaneous  equations  are  obtained 
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0 


■oox  i  +  a)  vsin6 
q  q  o 


( R^+R)  i  ^  +  a)  VC0S6 

oK  e  q  o 


(A2-1.4) 


which  in  turn  gives  the  following  nonlinear  differential  equations 


Pi 


fd 


( x  ,+x^ ) i ^  ,  x  ,,(R  +R)  x^(x  +0 

fdv  d  e'  fd  +  afdx  e  .  _  afdv  q  e'  ci+n)-j 


x  r  J  v 

+  — r—  sin6  -  E,  . 

A  ,  A,  fd 


rfd(xd+xe) 


(A2-1.5) 


rfd  xafd  .  .  xffd(Re+R)  .  xffd(xq+xe)  ,, 

P1d=— A—  ’fd  +  - ^ - ’d  -  H-— (l+n)iq 


X-P  r  i  V  r  n  I 

A  ,  A ,  fd 

d  d 
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^  0+n)ifd  -  d+n)id 
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Aq  q 


COS  6 


(A2-1.7) 


where 


xafd  ‘  xffd^xd+xe* 
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x  +xQ 

=  q 


00, 


Substituting  pi^  and  pi^  in  equation  (2.7)  gives 
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(A2-1.9) 


The  torque  and  voltage  regulator  equations  can  be  written  as 


"^in  xafd  „• 
pn  =  ^ —  — =? —  i  ^  i_  + 


(Xd-Xq>  .. 


T  T  'fd  'q  '  Tm  'dnq 
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(A2-1.10) 


P«  =  uon 


(A2-l.ll) 


PEfd  Tr  et  ’  Tr  Efd  +  ^Eo’Kr  etr^Tr  ”  us^ 


(A2-1.12) 


where 


2  2  .  2 
et  =ed+eq 


(A2-1 .12a) 


Equations  (A2-1.5)  to  (A2-1.7)  and  (A2-1.10),  (A2-l.ll)  can  be  grouped  as 
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APPENDIX  A2-2 


For  small  perturbations,  the  system  equations  can  be  lineari zed 
about  an  operating  point.  In  such  cases,  the  change  in  speed  deviation 
is  negligibly  small,  so  that  oo/a)Q  =  1.  To  include  the  voltage  regulator 
action  in  the  general  linearized  model,  it  is  necessary  to  obtain  an 
expression  for  perturbation  of  terminal  voltage  in  terms  of  the  states 
chosen. 


From  equation  ( A2- 1 ,12a) 


/ 


(A2-2.1) 


The  voltage  regulator  equation  in  the  linearized  form  is 


(A2-2.2) 


Substituting  the  linearized  forms  of  Ae^  and  Ae^  in  (A2-2.1)  one  gets 


(A2-2.3) 
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The  linearized  voltage  current  torque  relations  are  written  as 
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(A2-2.5) 

where  c. .  =  a--;  1  =  1.  2,  3 

*  J  '  U 

j  =  1,  2,  3 
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CHAPTER  3 


CLOSED  LOOP  QUASI-OPTIMAL  EXCITATION  CONTROL  FOR  A  SINGLE  MACHINE 

3.1  Introduction 

The  disturbances  appearing  on  a  power  system  are,  in  general, 
not  known  in  advance.  To  the  best  knowledge  of  the  author,  the 

available  methods  of  optimization  involve  iteration  even  for  a  known 

\ 

disturbance.  This  is  costly  in  terms  of  computing  time.  The  time 
involved  in  finding  the  optimal  control  by  digital  computers,  after  a 
fault  appears  on  a  system,  may  be  so  long  that  the  system  becomes  unstable 
before  a  control  strategy  is  determined.  A  suitable  control  must  overcome 
or  avoid  this.  Secondly,  a  predetermined  strategy  (found  by  iterative 
schemes)  will  not  be  optimal  for  other  than  the  design  disturbance.  A 
closed  loop  control  in  which  the  excitation  depends  on  actual  conditions 
and  requires  a  minimum  of  computation  in  its  implementation  is  desirable. 
This  chapter  discusses  such  a  control. 

The  stability  problem  for  a  single  machine  infinite  bus  system 
is  formulated  in  Section  2.  A  quasi -optimal  state  feedback  excitation 
control  for  the  linearized  system  is  suggested  in  Section  3.  Section  4 
investigates  the  closeness  of  the  quasi -optimal  solution  to  the  optimal 
one.  In  Section  5,  the  method  is  extended  to  the  nonlinear  machine  model. 
For  implementing  the  quasi -optimal  closed  loop  control  on  real  systems, 
a  practical  compromise  is  suggested  in  Section  6.  The  properties  of  a 
functional  which  depends  on  the  states  of  the  system  is  discussed  in 
Section  7.  The  results  obtained  for  the  single  machine  infinite  bus  system 


* 
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are  discussed  in  Section  8. 


3.2  Statement  of  the  problem 


Fig.  3.1  System  block  diagram 

Assuming  that  the  synchronous  machine  is  equipped  with  a  static 
excitation  system,  the  regulator  time  constant  is  very  small  with  respect 
to  the  machine  time  constants  and  initially  x^  (equation  2.19)  is  taken 
as  zero.  Then  the  voltage  regulator  equation  becomes  a  static  equation. 

The  exciter  in  such  a  case  can  be  driven  to  the  ceilings  almost  instant¬ 
aneously.  The  dynamic  equations  for  a  single  machine  case,  as  given  in  the 
previous  chapter,  are 

l(t)  =  f ( X.)  +  d  u(t)  (3.1) 

The  magnitude  of  the  field  voltage  is  bounded  by  the  exciter  ceiling  so 
that  the  normalized  scalar  control  u(t)  is  constrained  in  magnitude  by 
the  relation  |u(t)|<  1.  It  is  assumed  that  u(t)  is  a  piecewise  continuous 


function  of  time. 


I  *■  I 
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The  system  under  consideration  returns  to  the  stable  equil¬ 
ibrium  point  (see  the  definition  given  in  Chapter  1)  if  both  velocity 
and  acceleration  (n,  h)  of  the  machine  decreases  to  zero  following  a 
disturbance  while  the  final  rotor  angle  remains  constrained  in  the  range 
from  0  -  tt/2  radians. 

The  optimization  problem  can  be  stated  as:  Find  the  admissible 
control  u(t)  which  transfers  system  (3.1)  from  the  set  of  given  initial 
states  X(0)  to  the  desired  final  states 


n(tfT*4(tfT  =  0 


(3.2) 


So  that  the  cost  functional 


(3.3) 


o 


is  minimized,  since  it  is  desired  that  the  transients  are  eliminated  in 
the  least  possible  time. 

3.3  The  linearized  system 

The  linear  system  is  valid  for  only  small  disturbances.  For 
such  cases,  the  rotor  angle  does  not  traverse  far  beyond  the  operating 
point.  It  is  likely  that  the  rotor  angle  will  always  remain  between  0  and 
tt/2  radians.  For  stable  operation  following  a  disturbance,  it  is  desired 
that  the  machine  velocity  and  acceleration  should  drop  to  zero  in  the 


: 


smallest  possible  time.  Rewriting  the  linearized  machine  equations 
from  equation  (2.20)  (disregarding  the  regulator  time  constant)  leads 
to 
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AX  =  C  AX  +  d  u ( t ) 


(3.4) 


and  the  swing  equations  become 


A6  =  0J  n 

o 


(3.5) 


n  '  C41  A1fd  +  c42  A1d  +  c43  A1q 


(3.6) 


Step  1 


Differentiate  equation  (3.6)  with  respect  to  time  to  get 


n  =  c4lPAifd  +  c42  pAid  +  c43  PAiq 


(3.7) 


Substitute  the  equations  for  pAifd>  pAirf  and  pAiq  from  equation  (3.4)  into 
(3.7)  to  obtain 

n  =  (c41  C11  +  c42  C21  +  c43  C31 ^A1fd  +  *c41  C12  + 

c42  c22  +  c43  c32*A1d  +  *C41  C13  +  c42  c23  +  C43  c33^A1q 
+  ^c41  c15  +  c42  c25  +  c43  C35^A6  +  ^c41dl  +  c42d2^u^^ 


(3.8) 


or 


ri  =  L(Aifd,  Ald,  Ai  ,  AS)  +  b  u(t) 


(3.9) 


34 


Since  Aifd,  Aid,  Aiq  and  A6  are  functions  of  time,  the  above  equation  can 
be  written  as 


n  =  L(t)  +  b  u (t) 


(3.10) 


where 


L(t)  =  (c41  C]1  +  c42  c21  +  c43  c31 )Aifd  +  (c41  c12  +  c42  c22 


c43  c32)Aid  +  (c41  c13  +  c42  c23  +  c43  c33)a1, 
+  ^c41  C1 5  +  c42  c25  +  c43  'W4'5 


b  C41  dl  +  c42  d2 


(3.11) 


Equation  (3.10)  can  be  represented  by  the  block  diagram  in  figure  3.2. 


Fig.  3.2  Double  integral  plant  with  disturbance  L(t). 

This  is  a  so  called  double  integral  plant  where  the  input  to  the  plant  is 
modified  by  an  additive  term  L(t),  the  value  of  which  is  generally 
unknown.  (Properties  of  L(t)  are  given  in  Section  (3.7)).  The  system 
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(3.10)  is  called  controllable^  only  if 

| L ( t ) |  <  | b  u ( t) |  (3.12) 

or  | L ( t )/b |  <  1  (3.13) 

For  small  disturbances,  the  quantity  L(t)  is  small  enough.  It  is  assumed 
that  L(t)  is  at  least  piecewise  continuous. 


Step  2 


Fig.  3.3  A  portion  of  L(t)  and  control  u(t) 
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L(t)  depends  on  the  states  of  the  system  (3.4)  and  its  initial  value  is 
known.  Suppose  at  t  =  tQ,  L(t)  is  l_Q.  The  optimal  control  problem  can 
be  reformulated  as  follows: 

Given  the  system 


where 


n  =  %(t). 


u  •  < 

ormn  - 


u0(t)  =  L0  +  b  u(t) 


u0(t)  <  U 


omax 


(3.14) 

(3.15) 


Find  the  admissible  control  that  forces  the  system  (3.14)  from  any  initial 
state  {n/0j,  n^}  to  the  origin  in  shortest  possible  time. 

The  analysis,  using  Pontryagin's  minimum  principle  (P.M.P)  for 

T2l 

time  optimal  control  of  the  system  (3.14)  is  given  in  several  references1  J 
and  is  briefly  outlined  here  for  completeness  only. 

For  convenience,  assume  x-j  =  n;  then  system  (3.10)  can  be  rewritten 
as 


x-j  =  (3.16) 

*2  =  Uo(t) 

[These  x's  may  not  have  any  relation  with  X_  in  equation  (3.4)]. 

The  Hamiltonian  for  (3.16)  is  given  by 

H  =  1  +  x2(t)  p-j(t)  +  uQ(t)  p2(t)  (3.17) 

If  x  (t),  p  (t)  and  u  (t)  correspond  to  the  optimal  state,  costate  and 
control  respectively,  then  by  the  minimum  principle 

1  +  x*(t)  p*(t)  +  U*(t)  p^(t)  <  1  +  x*(t)  p*(t)  +  uo(t)  p*(t)  (3.18) 


37 


holds  for  all  admissible  uQ(t)  and  for  te[tQ,  tf].  Equation  (3.18) 
gives 


* 

u0(t)  ■ 


uomax>  if  p2(t)  "  0 


uomin’  if  P2(t)  "  0 


The  costate  variables  satisfy  the  relation 


A  t + \ - _  3H 

ax^t) 


=  0 


3H 


P2(t)  '  3x"(t)  '  Pl^ 


Let  it i  and  tt 2  be  the  initial  values  of  the  costates,  then 


p-j  (t)  =  tt-j  =  constant 


P2 ( ^ )  ^2  ”  ^1^ 


(3.19) 


(3.20) 

(3.21) 


(3.22) 

(3.23) 


Note  that  p2(t)  is  a  straight  line  in  the  p2  -  t  plane  and  can  intersect 
the  time  axis  once  only.  So  from  relation  (3.19),  the  optimal  control 
is  piecewise  constant  and  can  switch  at  most  once.  Assume  uQ(t)  =  a 
(which  is  either  uQmax  or  uomin)  and  let  x-j(o)  =  ^  and  x2(o)  =  e2- 
Equation  (3.16)  can  be  solved  with  these  initial  conditions  for  constant 
control  to  obtain  the  relations 
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x2(t)  =  +  at 

x -j  ( t )  =  s1  +  e2t  +  \  at2 

Next,  eliminate  t  to  find 

x  -  r  +  1_  x2  1_  P2 

X1  *1  2a  X2  2a  ^2 

where  t  =  (x2  -  ?2)/a 


(3.24) 

(3.25) 

(3.26) 
(3.26a) 


Equation  (3.26)  is  the  equation  of  the  trajectory  in  the  x-j  -  x2  plane 
originating  at  (?^,?2)  as  affected  by  the  action  of  the  control  a.  The 
switch  curve  is  the  locus  of  all  points  (x-j  ,x2)  which  can  be  forced  to 
zero  by  the  control  uQ  =  a  and  is  given  by 


If 


Y°  |(xi*x2)  :  xi  -  2[LQ  -  | b | Sgn{x2J J  0 

x2 

^  '  X1  [LQ  -  S b | Sgn{x2>] 


(3.27) 

(3.28) 


then  the  control  law  can  be  given  as 


* 

u0(t)  = 


u  •  if  y°  >  0 
oirnn  L 


u  if  y°  <  0 

omax  L 


(3.29) 


If  =  0,  then  u*(t)  is  uomin  for  x£  >  0,  otherwise  u*(t)  is  uomax.  But 

for  a  physical  system,  the  trajectory  does  not  continue  along  the  curve 
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1°  =  0  for  a  finite  interval  of  time^-L 
Step  3 

Given  the  system  of  equations  (3.16),  for  any  set  of  initial 
states  the  optimal  control  uQ ( t )  is  obtained  by  control  law 

(3.29).  The  value  of  the  optimal  control  u  (t)  is  decided  by  the  relation 

(3.15).  If  b  is  positive,  uQmax  implies  u  =  1  and  if  b  is  negative 

uomax  lmP^es  u  =  "1  and  vice  versa. 

The  term  L(t)  depends  on  Ai^,  Ai^ ,  Ai^  and  A6  (and  derivative 
of  the  disturbance  function,  if  there  is  a  sustained  disturbance).  The 
initial  value  of  L(t)  is  calculated  depending  on  the  initial  values  of 
currents  and  rotor  angle.  The  optimal  control  is  then  obtained  for  the 
constant  value  Lq.  This  control  u(t)  (*  for  the  optimal  values  are 
dropped  for  the  rest  of  the  analysis)  is  then  used  to  solve  the  system 
of  equations  (3.4)  for  a  small  time  step. 

Step  4 

At  the  end  of  the  first  interval,  at  t  =  t-j  (fig.  3.3),  suppose 

that  the  states  in  the  x-j  -  x2  plane  are  and  the  corresponding  values 

of  other  states  are  { Ai^^ ( i ) ,  A 1  d ( 1 )  ’  Alq(l)  anc*  A<5(1)^‘  Usin9  these 
new  values  of  currents  and  angle  recalculate  L(t).  Suppose  L(t)  =  L-| 
at  t  =  t-j.  If  L-j  =  L  ,  the  control  law  (3.29)  still  holds  good  and 
the  process  is  continued.  However,  if  L-|  f-  Lq,  then  the  second  problem 
is: 


Given  the  system  of  equations 


'• 

* 


i 
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Fig.  3.4  Moving  switch  curves  in  phase  plane. 
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n  =  L,  +  bu(t) 


=  u-,  ( t) , 


(3.30) 


Find  the  admissible  control  that  forces  the  system  (3.30)  from  the  given 
initial  states  to  the  origin  in  minimum  possible  time.  Proceeding 
as  before,  the  switch  curve  is  given  by 


(3.31) 


and  the  switching  function  and  control  scheme  is  similar  to  (3.28)  and 


(3.29)  respectively  with  LQ  replaced  by  L] .  Similarly  at  t=t  ,  the  switch 
curve  for  L(t)  =  L  is  given  by 

r 


Figure  3.4  shows  the  set  of  moving  switch  curves  in  the  phase  plane. 

3.4  Optimality  of  the  control 

It  can  be  seen  that  for  finding  the  control  u(t),  the  quantity 
L(t)  has  been  approximated  by  a  staircase  function.  Since  L(t)  is  being 
approximated,  it  would  appear  the  control  scheme  (3.29)  would  give  rise 
to  a  sub-optimal  solution.  In  fact,  if  equation  (3.10)  was  solved  by 
applying  this  control  scheme,  only  a  sub-optimal  solution  would  be  found, 
the  nearness  to  optimality  being  decided  by  the  size  of  the  step. 


However,  a  close  look  would  reveal  that  equation  (3.10)  is  just  a 


* 


■ 

- 
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dummy  equation  which  is  needed  to  find  the  optimal  control.  The 
information  obtained  from  (3.10)  is  deciding  the  control  at  the  sample 
points  { t i ,  — .ts>  ....}  which  in  turn  is  being  fed  to  the  system 

(3.4) .  It  will  be  shown  that  irrespective  of  the  variation  of  L(t)  in 
the  range  tc[t  ,  t  the  control  u  is  +1  (figure  3.3). 

Observe  that  at  t  =  t  ,  the  control  is  +1  in  figure  3.3  and  it 
is  optimal  since  by  proposition  L(t)  does  not  change  instantaneously. 

Again  at  t  =  t-j ,  the  control  is  +1 ; by  the  control  law  (3.29)  at  that 
instant  also  the  control  is  optimal.  Since  right  hand  side  of  equation 
(3.1)  is  linear  in  the  control  u,  the  time  optimal  control  for  system 

(3.4)  is  bang  bang  i.e.,  either  it  is  +1  or  -1.  The  only  other 
possibility  in  this  particular  case  is  that  the  control  switches  from  +1 
to  -1  and  then  switches  back  to  +1  in  the  interval  { t  , t -| } .  This  means 
that  the  control  switches  twice  in  the  first  time  step.  But  this  interval 
can  be  chosen  small  enough  so  that  the  control  really  cannot  have  two 
intermediate  switches.  Moreover,  the  number  of  switchings  for  system  (3.4) 
is  finite.  It  is  improbable  that  the  control  will  switch  twice  in  each 
subinterval.  It  can  be  seen  in  the  figure  that  in  the  intervals  (t  »t-j}, 

{t-,,t0},  . {t  9,t  -,},  irrespective  of  the  variation  of  L ( t )  the 

control  is  +1  i.e.,  in  this  interval,  the  staircase  approximation  of  L(t) 
does  not  affect  the  optimal  scheme. 

However,  the  control  switches  at  t  =  t  and  this  switching  is 
determined  by  a  constant  value  L$.  But  since  L(t)  does  not  remain  constant 
in  (t  ^t^},  the  correct  switching  can  be  anywhere  in  this  interval. 

The  error  in  switching  will  be  negligible  since  the  subintervals  can  be 
made  very  small  as  allowed  by  numerical  integration  procedures  with 


<  . 

> 
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sufficient  accuracy. 

The  control  obtained  by  the  above  scheme  is  optimal  over  each 
subinterval,  except  possibly  in  the  neighbourhood  of  the  switch  points. 

As  the  solution  is  time  optimal  in  each  subinterval,  it  is  optimal  over 
the  whole  interval  (A  similar  work  done  by  Oldenburger^  is  given 

in  Section  A3-1  of  appendix  A3). 

3.5  The  nonlinear  model 

For  large  disturbances,  the  linearized  machine  model  gives  a  poor 
approximation.  In  such  cases,  the  rotor  angle  might  swing  past  90  degrees 
and  there  is  a  possibility  that  the  machine  will  settle  momentarily  at 
the  unstable  equilibrium  point.  For  stable  operation  of  the  machine,  the 
velocity  and  acceleration  of  the  machine  should  drop  to  zero  following 
a  disturbance  while  the  final  rotor  angle  should  be  constrained  between 
0  and  tt/2  radians. 

In  finding  the  control  function,  it  is  recognized  that  n  is  very 
small.  So  the  terms  where  n  (speed  deviation)  is  coupled  to  other  system 
states  are  neglected.  (This  is  not  necessary  but  simplifies  the  algorithm). 
Differentiating  the  nonlinear  swing  equation  (A2-1.13)  yields 

P3(<5/«o)=  (a41  1fd  +  a42  id)Piq  +  a41  1q  pifd 

+  a42  ^q  Pn'd  (3.33) 

Substituting  the  nonlinear  voltage  current  relations  from  equation  (A2-1.13) 
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in  equation  (3.33)  gives 


a41  a31 

1  fd 

+  a42  a 

32  nd  + 

(a4l  a13 

+  a42  a23)iq 

+  (a41 

a32  + 

a42  a31 

^fd  jd + 

(a4i  an 

+  a42  a21  +  a41  a33) 

ifd  \ 

+  {a41 

a12  + 

a42  a22  + 

a42  a33^ 

•  • 

’d  ’q 

+  (a41 

nfd  + 

a42  nd^ 

aoc  cos 6 
35 

+  (a41  a 

15  +  a42  *25^  sin5 

+  (a41 

d-j  +  a 

'42  ^2^ 

q  Efd 

(3.34) 

Equation  (3.34)  can  be  written  as 


p'j(6/a)0)  =  L(t)  +  b ( t )  u (t) 


(3.35) 


2  2  2 
where  L(t)  =  a41  a31  ifd  +  a42  a32  id  +  (a41  a]3  +  a42  a23)iq 


+  (a41  a32  +  a42  a31)ifd  id  +  (a41  an  +  a42  a21  +  a41  a33) 


ifd  ’q  +  (a41  a12  +  a42  a22  +  a42  a33^d  ]q 


+  (a41  Vd  +  a42  V  a3  5COSS  +  *a41  al  5  +  a42  a2  5^  ’q  Sln6 


b(t)  -  (a4i  d-|  +  a42  d2)iq 


(3.36) 


« 
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Assume  that  L(t)  and  b(t)  are  constant  at  their  initial  values  LQ  and 
bQ  respectively  and  let  x-|  =  6/(joq,  then  (3.35)  can  be  written  as 


x2  x3 
x3  =  uo(t) 


(3.37) 


where 

and 


=  Lo  +  bo  “W 


u  .  <  u  (t)  <  u 

omin  -  ov  '  -  omax 


(3.38) 


The  optimization  problem  is  to  find  the  admissible  control  that 
forces  the  system  (3.37)  from  any  set  of  initial  states  U-,  (x-j  (o) ) , 

^ 2 ( x2 ^ 0 ) )  *  ^3^3(0))}  to 


0  <  x1(tf)  <  Tt/2co0 

x2(tf )  =  0  <3-39> 

X3(tf )  =  0 

in  the  shortest  possible  time. 

As  in  the  case  for  the  second  order  system  (3.16),  the  Hamiltonian 
can  beconstructed  and i t  can  be  shown  that  the  H-minimal  control  is  piecewise 
constant  and  can  switch,  at  most,  twice.  The  system  of  equations  (3.37) 
can  be  solved  for  constant  uQ ( t )  =  a  to  give 


Xg  =  at  +  £3 


(3.40) 
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'  E  Qt  t  I  *■  J.  I  _ 

v2  2  ?3r  ?2 


at3  ^3^ 

xlB“T-+-T-+52t  +  *l 


from  (3.40) 


t  =  (x3  -  53)/o 


Substituting  for  t  in  (3.41)  one  gets 


x0  =  Co  + 


x  2  C  2 
x3  ^3 


2  ^2  2a  2a 


X1  = 


a 

[X3  '  ?3\3  ,  ^3 

(x3  -  hY  + 

x3  -  s3) 

6 

|  2 

a  |  52 

a 

+  5. 


This  leads  to 


X  =  J_£  3  .  f2fi  +  +  J_x  3  ^3  (  _  ]_  x  2 

°  a  51  ,  2  X3  a  V?2  2a  X3 


1  .  2  ^3 

oa 


6a 


(3.41) 

(3.42) 

(3.43) 


(3.44) 

(3.45) 


(3.46) 


Substituting  (3.44)  in  this  expression  yields 


Or 


1 


3  5253 


1  .2-3 

oa 


a 


+  5 


1  3a2  3 


1  +  X2X3 


a 


^  1  3  X2X3  .  „ 

X1  .  2  X3  ‘  a  C1  .  2  -3 


1  .  3  5253 


3a 


3a 


a 


(3.47) 


(3.48) 


Expression  (3.48)  is  the  equation  of  a  surface  in  the  x^-x^-x^  plane  for 
any  initial  state  (^*|  »^2 *^3 ^ ’  T*ie  ecluatl'on  °f  the  trajectory  passing 
through  x-j  =  tt/2u)o,  x2  =  0  and  x^  =  0  is  given  by 


•  • 
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x0x0 
it  _  2  3 

2co_  a 

0 


1 


(3.49) 


which  is  one  of  the  switching  surfaces  in  the  x-j-x^-x^  plane. 

The  other  surface  which  passes  through  the  origin  is  given  by 


where 


X2X3 

a 


a  =  L_  +  b_  Sgn 
o  o 


x 


2  ‘  2[Lq  +  bQSgn{x3}]  f*  o 


0 


(3.50) 

(3.51) 


In  the  region  between  these  two  surfaces,  the  control  is  decided  by  the 
states  x^-x^  and  the  switching  curve  is  given  by 


Xo  - 


x. 


2  L  L0  +  boSgn(x3}J 


=  0 


b„  <  0 
o 


(3.52) 


The  following  steps  are  involved  in  finding  the  control 


1.  Determine  x^  and  x^. 

2 

x3 

2.  Determine  I  =  x2  -  2[l  +  b  Sgn{x,}J  ;  bo  <  0  (3'53) 

L  0  0  3 

3-  If  l  >  °>  “  =  %min  otherwise  a  =  uomax' 

4.  If  l  -  0  and  x3  >  0,  a  =  uQmin  otherwise  a  =  uomax. 

5.  Determine  V  =  x,  -  ^ - f1  +  “T  V*  (3’54) 

I  0  3a 

y  -  x  _  X2  3  +  ,  x  3  /g 

h  "  X1  a  3a2  X3 

If  J  <0  and  l  >  0 
1  '  2 


6. 


. 


u(t)  =  Sgn 


(3.56) 


x 


3 _ 

2  "  2 L Lq  +  boSgn{x3}j^ 


b  <  0 
o 


7.  If  l  >  0,  u(t)  =  1  ;  b  <  0 

1 

If  £  <0,  u ( t )  =  -1  ;  b  <  0 
2  0 


(3.57) 


The  system  does  not  continue  along  the  switching  surfaces  for  a  finite 

m 

interval  of  time  due  to  system  imperfections'-  J.  The  process  outlined 
above  is  continued  until  the  desired  final  states  are  reached. 

3.6  Proportional  Control 

So  far  it  has  been  considered  that  the  regulator  time  constant  xr 

is  zero,  so  that  the  exciter  can  be  driven  to  the  ceiling  voltages 

instantaneously.  This  means  that  for  small  deviations  of  velocity  and 

acceleration,  the  exciter  will  still  be  repeatedly  driven  to  the  limits. 

This  causes  unnecessary  wear  and  tear  to  the  equipment  and  is  equivalent 

to  'firing  a  cannon  to  kill  a  bug1.  It  deteriorates  the  machine  internal 

voltage  and  worsens  the  terminal  voltage  of  the  machine.  Moreover,  the 

reproduction  of  complicated  control  functionsby  physical  equipment  in 

[4] 

industry  is  generally  costly  and  practical  compromises  must  be  made  . 

To  allow  normal  voltage  regulator  action,  a  small  deadzone  has  to  be 
provided  so  that  when  the  oscillations  are  small  enough,  the  voltage 
regulator  takes  over.  It  is  usually  not  convenient  in  practice  to  use 


u(t)  =  Sgn  l 


(3.58) 


as  the  control  equation  but  to  employ  the  proportion  relationship 
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u  =  k  l 


(3.59) 


where  k  is  as  large  as  possible.  For  k  =  «,  equation  (3.59)  reduces  to 
(3.58).  An  optimum  value  of  k  has  to  be  found  by  measurement  or  computer 
simulation.  Too  high  a  gain  throws  the  system  into  a  steady  hunting. 
Equation  (3.59)  gives  approximate  bang  bang  control  for  large  disturbances. 
When  the  oscillations  gradually  die  out,  normal  voltage  regulator  action 
takes  over. 

For  small  disturbances  the  acceleration  term  in  l  is  negligible 
compared  to  the  velocity  term,  so  that  proportional  control  gives 


u(t)  =  k l 


(3.60) 


=  kn 


which  is  the  velocity  signal  as  used  at  present.  For  large  disturbances 
a  combination  of  velocity  and  acceleration  terms  is  determined  by  L(t), 
i.e.,  by  the  various  currents  and  the  rotor  angle.  References  [5, 6, 7, 8] 
bear  this  out.  The  proportional  control  is  not  quasi -optimal  and,  in 
general,  quasi-optimal  or  optimal  control  may  not  be  impl ementable. 

3.7  Some  comments  on  L(t) 

Rewriting  equation  (3.9) 


h  =  L(Aifd,Aid,Aiq,A6)  +  b  u(t) 


(3.61) 


The  arguments  of  L  are  functions  of  time,  so  without  loss  of  generality 


;  6U0  ‘  '  ir<1  1  *  ^  **  1  *  j 
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equation  (3.61)  may  be  expressed  as 

n  =  L ( t )  +  b  u(t),  | u ( t) |  <  1  (3.62) 

The  block  diagram  for  system  (3.62)  is  given  in  figure  3.2. 

The  concept  of  controllability,  as  given  in  reference  [2]  is  that 
a  system  is  controllable  if  it  is  possible  to  drive  any  state  of  the 
system  to  the  origin  in  finite  time.  Here  controllability  means  complete 
controllability  (global  controllability).  For  example,  the  system 

n  =  u(t) ,  | u ( t ) |  <  1  (3.63) 

is  completely  controllable.  But  if  the  upper  and  lower  bounds  of  the 
control  u(t)  is  such  that  0  <  u(t)  <  1,  then  system  (3.63)  is  not  completely 
controllable  since  there  are  some  states  in  the  state  space  which  can  never 
be  brought  to  the  origin  with  the  application  of  this  control.  It  is 
obvious  that  the  system  (3.62)  is  controllable  only  if^ 

| L ( t ) / b |  <  1  (3.64) 

The  term  L(t)  could  be  considered  an  additive  disturbance  term  to 
the  input  of  system  (3.63)  which  is  governed  by  a  set  of  differential 
equations.  It  can  be  seen  that  L(t)  does  not  contain  any  output  or 
derivative  of  the  output  feedback  term  and  also  L(t)  can  have  any  value 
satisfying  constraint  (3.64). 

For  the  nonlinear  system  (3.35)  L(t)  contains  feedback  terms  such  as 


- 
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6  and  n.  Moreover,  for  large  disturbances,  it  may  not  satisfy  the 
constraint  (3.64)  for  all  t^:[to,t^].  In  such  cases,  the  solution 
obtained  with  the  control  scheme  is  not  optimal. 

To  summarize,  L(t)  should  have  the  following  properties,  so  that 
the  solution  will  be  closest  to  optimal. 

1.  It  should  be  continuous  in  the  time  interval,  at  least  it 
should  not  have  any  jump  discontinuities. 

2.  It  should  satisfy  the  relation 

|  L(t)/b |  <  | u ( t )  |  for  all  te[tQ,tf]  (3.65) 

3.  Output  feedback  terms  should  not  dominate  L(t)  and  L(t) 
should  not  have  any  other  constraint  except  (3.65). 

A  sub-optimal  solution  is  obtained  if 

1.  L(t)  contains  output  or  derivatives  of  output  feedback 
which  are  dominant. 

2.  The  uncontrollable  region  is  followed  by  a  long  controllable 

-  [4] 

section1-  . 

An  example,  showing  that  a  sub-optimal  solution  is  obtained  when 
L(t)  contains  output  terms,  is  given  in  Section  A3-2  of  appendix  A3. 

3.8  Discussion  of  results 

A  power  system  comprised  of  synchronous  generators  and  transmission 
network,  whose  parameters  are  given  in  section  A3-3  of  appendix  A3,  was 
simulated  on  a  digital  computer.  Damper  windings  and  saturation  of  the 
machine  and  governor  action  are  not  considered.  The  system  configuration 
is  given  in  figure  2.2. 
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i)  The  linearized  model  -  10%  torque  step. 

For  small  disturbances,  the  nonlinear  differential  equations  of 
the  system  can  be  linearized  about  an  operating  point.  A  10%  input 
torque  step  is  considered  for  such  a  system. 

Figure  3.5  shows  the  effect  of  exciter  ceiling  on  the  response 
of  the  system.  With  the  bang  bang  control  the  time  for  stabilization  and 
the  angular  deviation  increases,  as  the  exciter  ceiling  is  decreased; 
that  is,  the  effectiveness  of  the  control  decreases. 


Fig.  3.5  Rotor  angle  time  characteristics  for  10%  torque  step 

(Bang  bang  control). 

1.  Ceiling  voltage  ±3  p.u. 

2.  Ceiling  voltage  ±5  p.u. 


* 
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The  optimal  bang  bang  control  results  in  a  deadbeat  response. 

While  the  use  of  proportional  control  provides  a  well  damped  stabilization, 
it  is  not  optimal  and  takes  more  time  to  stabilize.  Voltage  regulator  action 
alone  leads  to  dynamic  instability.  These  are  shown  in  figures  3.6,  3.7 
and  3.8. 


Fig.  3.6  Angle  time  characteristics  for  10%  torque  step. 

1.  Optimal  control. 

2.  Proportional  control. 

3.  Voltage  regulator  alone. 


ACCELERATION 


54 


-0.001 


Fig.  3.7 


Velocity  vs  acceleration  plot  corresponding  to  Fig.  3.6. 

1.  Optimal  control. 

2.  Proportional  control. 
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Fig.  3.8  Angle  vs  velocity  plot  corresponding  to  Fig.  3.6. 

1.  Optimal  control. 

2.  Proportional  control. 


One  factor  of  interest  in  stability  studies  is  the  gain  of  the 
voltage  regulator.  The  regulator  gain  was  varied  from  -4.5  to  -1000 
with  the  results  shown  in  Table  3.1.  It  is  found  that  a  gain  of  -1000 
drives  the  exciter  almost  in  bang  bang  fashion  even  for  a  small  disturbance. 
With  the  proportional  control,  a  higher  gain  regulator  takes  less  time  to 
stabilize.  A  gain  of  -100  appears  to  be  optimal  in  that  the  exciter  is  not 
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driven  too  hard  but  effective  stabilization  is  obtained.  Increasing  the 
regulator  gain  implies  a  shift  of  the  roots  of  the  linearized  system  to 
the  right  hand  side  of  the  complex  frequency  plane.  The  bang  bang  control 
assumes  that  the  regulator  gain  is  infinite. 


Case 

Reg .Gain 

Control  Signal 

Ceil ing 
Vol tage 
(p.u. ) 

Stabil  ity 

Time  to 
stabilize 
(sec) 

Final  rotor 
angle 

1 

-4.5 

Proportional 

±5 

Stabl e 

- 

68° 

2 

-4.5 

Voltage 

Regul ator 
only 

±5 

Unstabl e 

3 

-10 

Prop 

±5 

Stable 

- 

68° 

4 

-10 

V.R. 

±5 

Unstable 

- 

- 

5 

-100 

Prop 

±5 

Stable 

- 

68° 

6 

-100 

V.R. 

±5 

Unstable 

- 

- 

7 

-1000 

Prop 

±5 

Stable 

- 

68° 

8 

-1000 

V.R. 

±5 

Unstable 

- 

- 

9 

- 

Bang  bang 

±5 

Stable 

,18 

62.1° 

10 

— 

Bang  bang 

±3 

Stable 

.35 

65.7° 

Table 

3.1  10%  Torque  Step, 

1 inearized 

system. 

ii)  Nonlinear  model  -  30%  torque  step. 


This  is  a  relatively  large  disturbance  and  the  system  has  to  be 


r 

V 

;  ' 
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represented  by  non-linear  equations.  The  term  L(t)/b(t)  exceeds,  in 
magnitude,  the  control  u(t)  for  a  very  small  period  of  time.  According 
to  the  definition  of  controllability  introduced  in  section  3.7,  the  system 
is  uncontrollable  for  that  period  of  time. 


Fig.  3.9  Angle  time  characteristics  for  30%  torque  step. 

1 .  Bang  bang  control . 

2.  Proportional  control. 


900 


58 


Proportional  control. 
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Fig.  3.11  Angle  vs  velocity  plot  corresponding  to  Fig.  3.9. 

1.  Bang  bang  control. 

2.  Proportional  control. 


Fig.  3.12  Terminal  voltage  characteristics  for  30%  torque  step. 

1.  Bang  bang  control. 

2.  Proportional  control. 
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The  bang  bang  control  obtained  by  the  closed  loop  scheme  is  sub- 
optimal.  This  results  in  an  overshoot  in  the  angle  time  characteristics 
(fig.  3.9).  The  time  required  for  stabilization  is  fairly  small.  The 
proportional  control  leads  to  a  heavily  damped  oscillatory  response.  The 
corresponding  phase  plane  plots  are  shown  in  figures  3.10  and  3.11.  Figure 
3.12  shows  the  variation  of  terminal  voltage  with  time  for  the  bang  bang 
as  well  as  proportional  control.  Proportional  control  is  better  in  this 
regard. 

Hi)  Nonlinear  model  -  100%  torque  pulse  for  3  cycles. 

This  is  a  severe  disturbance.  The  term  L(t)/b(t)  remains  very 
large  with  respect  to  u(t)  for  a  considerable  period  of  time.  The  system 
initially  is  uncontrollable.  However,  this  uncontrollable  section  of  L(t) 
is  followed  by  a  long  controllable  section,  giving  a  sub-optimal  solution. 
The  angle  time  characteristics  (fig.  3.13)  show  a  large  overshoot  followed 
by  an  undershoot.  Proportional  control  produces  a  highly  damped  response 
with  a  longer  settling  time.  The  two  curves  are  nearer  to  each  other  than 
in  the  30%  torque  step  case.  In  figure  3.13  it  is  observed  that  the 
proportional  control  gives  a  lower  first  peak  than  the  bang  bang  control. 

The  bang  bang  control  causes  a  large  initial  slope  in  the  states  and  when 
the  disturbance  ceases  after  3  cycles,  leads  to  a  large  peak  in  angle  due 
to  inertia.  On  the  other  hand,  the  proportional  control  is  slower  and 
by  the  time  it  becomes  effective,  the  disturbance  is  not  in  effect, 
requiring  less  effort  and  consequently  the  first  peak  is  lower.  The  phase 
plane  plots  are  shown  in  figures  3.14  and  3.15.  If  in  this  case,  the  torque 


- 


ANGLE  ( deg.) 
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pulse  were  150/  for  3  cycles,  the  system  would  be  uncontrollable  for  the 
particular  ceiling  voltage. 


1 .  Bang  bang  control . 

2.  Proportional  control. 
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Fig.  3.14  Velocity  vs  acceleration  plot  corresponding  to  Fig.  3.13. 

1.  Bang  bang  control. 

2.  Proportional  control. 


SPEED  DEV.  (p.u.) 
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Fig.  3.15  Angle  vs  speed  deviation  plot  corresponding  to  Fig.  3.13. 

1 .  Bang  bang  control 

2.  Proportional  control. 


TERMINAL  VOLTAGE 
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Fig.  3.16  Terminal  voltage  characteristics  for  100%  torque  pulse. 

1.  Bang  bang  control. 

2.  Proportional  control. 


, 


FIELD  VOLTAGE  (p.u.) 
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Fig.  3.17  Field  voltage  variation  for  100%  torque  pulse. 

1 .  Bang  bang  control . 

2.  Proportional  control. 
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A  factor  of  interest  is  the  terminal  voltage  of  the  machine  (fig.  3.16). 

As  would  be  expected,  the  bang  bang  control  deteriorates  the  voltage 
regulation,  since  the  voltage  regulator  action  is  completely  overridden  by 
the  bang  bang  control.  This  is  offset  to  some  extent  by  faster  elimination 
of  the  transients.  Figure  3.17  shows  multiple  switching  in  the  field 
voltage.  This  is  partly  due  to  sudden  variations  in  L(t)  and  b(t)  and 
partly  because  of  quasilinearization  of  these  qualities.  This  is  dealt 
with  in  more  detail  in  the  next  chapter. 

Apart  from  30%  torque  step  and  100%  torque  pulse,  several  other 
cases  have  been  considered.  These  are  summarized  in  Table  3.2.  With  the 
voltage  regulator  alone,  all  cases  are  unstable. 


Case  Mode  of 

disturbance 

Control 

Signal  Stability 

Time  to 
stabil ize 
(sec) 

Final  rotor 
angle 

1 

20%  torque  step 

Bang  bang 

Stabl e 

.355  sec 

70° 

2 

20%  torque  step  Proportional 

Stable 

- 

74° 

3 

30%  torque  step  Bang  bang 

Stabl e 

.48 

84° 

4 

30%  torque  step  Proportional 

Stable 

- 

82° 

5 

50%  pulse  for 

3  cycles 

Bang  Dang 

Stabl e 

.668 

59.5° 

6 

50%  pulse  for 

3  cycles 

Pro portional 

Stabl e 

- 

60° 

7 

100%  pulse  for 

3  cycles 

Bang  bang 

Stabl e 

.988 

59.5° 

8 

100%  pulse  for 

3  cycles 

Proportional 

Stable 

- 

60° 

Table  3.2  Nonlinear  system,  torque  step  and  pulses 

(K  =  -  100,  celing  ±  5  p.u.) 


. 

. 


• 
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iv)  Nonlinear  model  -  Three  phase  fault  at  intermediate  bus;  cleared  in 
3  cycles  by  opening  a  line. 

The  bang  bang  control  gives  one  overshoot  and  a  small  undershoot 
(fig.  3.18).  While  the  fault  is  on  and  a  small  period  afterwards,  the 
system  is  uncontrollable.  If  the  final  rotor  angle  is  considered  free, 
the  system  attempts  to  stabilize  at  the  unstable  equilibrium  point  and 
then  slowly  creeps  further  out  of  phase.  The  constraint  on  the  final 
rotor  angle  prevents  this.  With  the  proportional  control,  the  system  has 
a  heavily  damped  stable  response  and  on  the  voltage  regulator  control  alone, 
is  unstable  (first  swing).  The  phase  plane  plot  (Fig.  3.19)  shows  that 
the  uncontrollable  period  continues  until  the  maximum  velocity  is  reached 
and  the  loops  produced  are  apparently  due  to  harmonic  torques.  Figures 
3.20  through  3.25  show  the  plots  of  velocity  vs  angle  phase  plane,  field 
and  terminal  voltage  characteristics,  field  and  different  armature  currents. 
Figure  3.26  shows  the  angle  time  characteri sties  for  a  ceiling  voltage 
of  ±  4.3  p.u.  A  ceiling  voltage  lower  than  this  causes  instability  for 
the  particular  disturbance. 

The  three  phase  short  circuit  is  also  considered  for  a  relatively 
lightly  loaded  condition.  This  and  the  summary  of  the  case  considered  is 
given  in  Table  3.3. 


*« 

» 


ANGLE  (deg.) 
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Fig.  3.18  Angle  time  characteristics  for  3  phase  fault. 

1 .  Bang  bang  control . 

2.  Proportional  control. 

3.  Voltage  regulator  only. 
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Fig.  3.19  Velocity  vs  acceleration  plot  corresponding  to  Fig.  3.18 


for  bang  bang  control  only. 


SPEED  DEV.  (p.u.) 
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Fig.  3.20  Angle  vs  velocity  plot  corresponding  to  Fig.  3.18. 

1.  Bang  bang  control. 

2.  Proportional  control. 


FIELD  VOLTAGE 
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Fig.  3.21  Field  voltage  characteristics,  proportional  control. 


FIELD  CURRENT 
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Fig.  3.22  Field  current  time  characteristics,  proportional 


control . 
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proportional  control . 
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Fig.  3.24  Quadrature  axis  armature  current  - 


3.2  4.0 


time  characteristics 


proportional  control. 


TERMINAL  VOLTAGE 
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Fig.  3.25  Terminal  voltage  variation  with  time,  proportional 


control . 
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Fig.  3.26  Angle  time  characteristic  for  the  minimum  allowable 


ceiling  voltage. 
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Operating 

point 

Case 

Control  Ceiling 

signal 

Stabi 1 i ty 

Time  to 

stabi 1 i ze 

Final  rotor 
angle 

(sec) 

o 

O 

II 

O 

«o 

1 

Bang  bang  ±7 

Stable 

1.011 

72° 

Pout  =  9235 

2 

Bang  bang  ±4.3 

Stable 

2.04 

o 

OJ 

CO 

3 

Proportion  ±5 
-al 

Stabi e 

- 

o 

LO 

CO 

4 

Voltage  regul¬ 
ator  only  ±5 

Unstable 

- 

- 

o> 

O 

II 

cn 

o 

o 

5 

Bang  bang  ±5 

Stable 

.694 

83° 

^out  73% 

6 

Proportion  ±5 
-al 

Stabi e 

— 

68° 

Table  3.3  Three  phase  short  circuit  cleared  in  3  cycles 
by  opening  a  1 ine 


< 
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APPENDIX  A3-1 


Oldenburger 


considered  a  system 


x  =  u(t)  -  L(t) 


(A3-1.1) 


where  L(t)  is  an  arbitrary  disturbance  term  and  the  control  is  constrained 
in  magnitude  by  | u | <  1 . 

By  a  geometrical  construction  and  also  applying  direct  method  of 
optimization,  he  found  a  switch  curve 


x  I  x  I 

Y  =  {(x-|,x2)  :  x1  +  2[-f- +\'Sgn{x2}]  =  0 


(A3-1 .2) 


Let  equation  (A3-1.2)  be  investigated  now.  For  L(t)  =  Lq ,  at  t  =  t 


uomax  ^  "  ^o 

u  .  =  -(1  +  L  ) 

omm  v  o' 


which  means  that  for  X2>0 


r° 

I  =  X 


1  2u 


and  for  X2<0,  l  =  x-j  - 


omin 
x  2 


2u 


omax 


(A3-1.3) 


(A3-1.4) 

(A3-1.5) 
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Equation  (A3-1.4)  and  (A3-1.5)  are  identical  to  equation  (3.45).  He 
considered  the  solution  to  system  ( A3 -1.1)  using  optimal  scheme  found 
from  equation  (A3-1.2)  to  be  suboptimal. 
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APPENDIX  A3-2 


This  section  considers  an  example  where  the  term  L(t),  considered 
in  Section  3.3,  contains  output  feedback.  The  problem  is  solved  by  the 
technique  suggested  in  Section  3.3  and  the  switch  curves  obtained  are 
compared  with  those  obtained  by  a  standard  optimization  method. 

The  equation  of  the  harmonic  oscillator  is 


..  9 

x  +  go  x  =  u(t) 


(A3-2.1) 


The  control  u(t)  is  constrained  in  magnitude  such  that  |u(t)|<  1.  The 
problem  is  to  find  the  admissible  control  which  transfers  system  (A3-2.1) 
from  the  given  initial  states  to  the  origin  in  smallest  possible  time. 

Considering  x  =  x-j  ,  x  =  x^  the  first  set  of  switch  curves  for  the 
system  (A3-2.1)  obtained  by  applying  Pontryagin's  Minimum  Principle 
(Reference  [2])  in  the  range  | x-j  | <  2  is 


(A3-2.2) 


(A3-2.3) 


which  are  two  semicircles  shown  in  figure  A3-2.1  (curve  1). 

2 

Next,  consider  -go  x  =  L(t),  then  equation  (A3-2.1)  can  be  written 


as 


X1  x2 

x2  =  L(t)  +  u(t) 


(A3-2.4) 


* 


. 
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Constraining  L(t)  in  magnitude,  so  that 

| L( t) |  <  1  (A3-2.5) 

or  |oo2x|  <  1  (A3-2.6) 

and  applying  the  technique  presented  in  Section  3.3,  the  switch  curve 
obtained  is 

x  J  xJ 

x-,  +  - -  =  0  (A3-2.7) 

2[1  +  w  x-|Sgn{x2)] 

For  x2  >  0  and  w  =  1,  the  equation  of  the  switch  curve  is 

X  2 

xi  +  2TT^T  =  0  (A3-2-8) 

or  Zx-f  +  2x1  +  x22  =  0  (A3-2.9) 

which  is  curve  2  in  figure  A3-2.1  in  the  range  |x-j|  <  1.  By  the  constraint 
(A3-2.6)  no  information  is  available  for  | x -j  |  >  1. 

Here  the  quantity  L(t),  i.e.  -x,  does  not  satisfy  the  restrictions 
imposed  on  it  in  Section  3.8.  (a)  Here,  L(t)  is  a  pure  output  feedback 

term,  but  as  indicated,  L(t)  should  not  be  dominated  by  output  feedback 
terms,  (b)  The  final  value  of  L(t)  is  zero  in  this  case,  which  puts 
an  extra  restriction  on  it  in  addition  to  (A3-2.5). 

Still  it  can  be  seen  from  figure  A3-2.1  that  for  states  sufficiently 
near  the  origin,  the  scheme  (A3-2.7)  gives  a  near  optimal  solution  implying 
that  if  L(t)  is  small  compared  to  u(t),  the  solution  will  be  close  to 
optimal.  It  can  be  seen  that  for  u>  <  .5,  the  switch  curves  (1)  and  (2)  in 
figure  A3-2.1  are  very  close  to  each  other  for  ] x-j  |  <  1.  Switch  curve  3 
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Fig.  A3-2.1  The  harmonic  oscillator  and  sub-optimal  control. 

F21 

gives  another  sub-optimal  solution1-  J  for  system  A3-2.1. 

From  the  above  analysis,  it  can  be  seen  that  as  long  as  the  term 
L(t)  is  constrained  in  magnitude,  the  control  scheme  gives  a  reasonably 
good  sub-optimal  solution  even  if  it  contains  output  feedback.  Of  course, 
sub-optimality  of  the  scheme  does  not  effect  the  solution  of  the  original 
power  system  equations  as  explained  in  Section  3.4. 


♦ 
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APPENDIX  A3-3 


Hamming's  modified  predictor  corrector  method  was  used  to  solve 
the  differential  equations  for  the  single  machine  case.  The  starting 
values  were  found  by  Runge-Kutta  Gill  numerical  integration  routine. 

The  data  for  the  machine  was  taken  from  Roy^^.  The  system  configuration 
is  given  in  figure  2.2. 


Machine  rating  46  MVA,  13.8  kv 


afd 

.562 

xtie  =  .7/section 

d 

.6803 

rtie  =  -07/section 

q 

.4363 

xtransf  " 

fd 

.000508 

00 

• 

II 

o» 

X 

ffd 

.692 

K  =  -07 

e 

3 

=  .0125  sec 

.00435 

Kr  =  -4.5  to  -1000 

The  operating  point  for  the  linear  model  (underexcited)  are 


Vdo 

ndo 


1.77935 

.30167 


=  60° 
o 

v  =  1 


.718637 


.662 
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et0  =  .85077  Qq  =  .0147 

Efd=  1 

The  operating  points  for  the  nonlinear  model  (normally  excited) 
are 


fdo 

=  2.473309 

6o  = 

60° 

do 

=  .564341 

V 

1 

qo 

=  .734433 

P  = 

.916 

0 

to 

=  1.05208 

II 

o 

O' 

.3324 

fd 

=  1.39 

CHAPTER  4 


EVALUATION  OF  A  QUASI-OPTIMAL  STATE  FEEDBACK  EXCITATION  SCHEME 

4.1  Introduction 

As  established  in  the  previous  chapter,  the  quasi -optimal 
control  is  optimal  for  the  linearized  system  model.  Use  of  this  model 
is  clearly  justified  when  L(t)  is  small  compared  to  the  control  u(t). 

For  large  disturbances  L(t)  is  quite  appreciable,  however,  and  the 
quasi -optimal  scheme  may  give  a  solution  quite  different  from  the  optimal 
ones.  In  order  to  evaluate  the  effectiveness  of  the  quasi -optimal 
scheme,  the  results  obtained  by  this  method  are  compared  with  those 
obtained  by  a  standard  optimization  technique. 

The  optimal  control  problem  stated  in  Section  (3.2)  of  Chapter  3 
is  solved  in  Section  2  using  the  steepest  descent  method.  The  results 
obtained  with  the  quasi -optimal  scheme  (closed  loop  scheme)  are  compared 
with  those  obtained  by  the  steepest  descent  method  in  Section  3.  Section 
4  investigates  the  effect  of  time  delays  on  the  optimal  and  quasi -optimal 
schemes . 

4.2  Numerical  example  using  steepest  descent 

n  pi 

The  steepest  descent  algorithm  constructed  by  Bryson  and  DenhamL  *  J 

T3l 

and  later  modified  by  VachinoL  J  to  take  into  consideration  control  problems 

★ 

with  inequality  constraints  has  been  used  here  .  For  starting  the  algorithm 
a  first  guess  for  the  optimal  control  (which  is  parametized  by  the  switching 

*  Details  of  the  steepest  descent  technique  is  given  in  Section  A4.1  of 
Appendix  A4. 


. 
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times  for  this  particular  problem)  and  final  time  is  required.  The 
control  provided  by  the  quasi -optimal  scheme  is  very  effective  in  this 
regard;  in  general,  it  produces  much  more  rapid  convergence  to  the  optimal 
solution  than  controls  only  slightly  different  from  it. 

i)  Linearized  model 

The  linearized  system  equations  rewritten  from  chapter  2 
(dropping  A  for  convenience)are 

X=CX+du(t)  (4.1) 

The  operating  values  of  u(t)  considered  for  linearization  are  1  per  unit. 
For  small  disturbances  the  exciter  may  not  be  driven  too  hard  and  it  is 
assumed  that  it  will  be  driven  to  ±  3  per  unit  only,  i.e., 

-4  <  u ( t )  <  2  (4.2) 

The  terminal  constraints  are  that  velocity  and  acceleration  of  the  machine 
at  t  =  tf  are  zero,  that  is, 

^1  =  "  x4  =  0  (4.3) 

^2  =  ”  x4  C41  x]  "  c42  x2  ~  c43  x3  _  ^  (4.4) 

Usually,  one  of  the  terminal  conditions  is  used  as  a  stopping 
condition  of  the  form  ft[x(tf),  tf]  =  0.  This  is  not  done  here  since  both 
velocity  and  acceleration  of  the  machine  are  alternating.  Satisfying 


+  Different  symbols  used  in  this  section  are  defined  in  appendix  A4. 
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only  one  of  the  terminal  conditions  does  not  necessarily  mean  that  the 
process  should  be  discontinued  there.  For  such  problems,  the  terminal 
time  tf  is  generally  chosen  as  a  stopping  condition^.  The  time 
required  for  stabilization  by  the  quasi -optimal  scheme  is  used  as  the 
stopping  time.  Assuming  that  the  quasi -optimal  scheme  only  gives  a  sub- 
optimal  solution,  the  optimum  time  will  be  less  than  that  obtained  by  it. 
The  stopping  condition  for  a  10%  torque  step  is  chosen  as 

n  =  .35  -  tf  (4.5) 

The  performance  index  is 

<f>  =  -  t-f  (4.6) 

Evaluating  the  values  of  the  parameters  q.  at  this  operating  point,  the 
terminal  conditions  are 

*1  =  x4(tf)  =  0  (4.7) 

i|>2  =  .06708  x1  (tf )  -  .02913  x2(tf)  +  .15437  x3(tf)  =  0 

(4.8) 


Assuming  two  switchings  of  the  control,  the  nominal  control  is 

u0(t)  =  2  t0  <  t  <  t, 
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uQ ( t)  =  -4  t-|  <  t  <  t2  (4.9) 

uQ(t)  =2  t2  <  t  <  tf 


Fig.  4.1  System  equations 


The  rewritten  system  equation  is 


X  =  f  i  [1 -h(t-t-| )]  +  f^C^lt-t-j )  -h(t-t2]  +  t^Chtt-^)  -  h(t-tp)] 


(4.10) 
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which  yields  the  equations  of  variation 

6X_  =  6f i [1  -  h(t-t-j)]  +  6f2[h(t-t  )  -  h ( t-t^) ] 

+  SfgChU^)  -  h(t-t^)]  +  (f  |-f  2)A(t-t-j  )6t-| 

+  (f2-f3)A(t-t2)8t2  +  f3A(t-tf)6tf  (4.11) 


This  vector  equation  can  be  written  in  the  component  form 


and 


6  X  =  F  SX  +  L6io 

X 

X  =  -FTx  X 


where  F  =  C 

A 


(4.12) 

(4.13) 

(4.14) 


L 


6d -j  A ( t-ti ) 
6d2A(t-tf) 
0 
0 
0 


-6d-|  A(t-t2) 
-6d2A ( t- t2 ) 
0 
0 
0 


2d  -j  A(t-tf ) 
2d2A(t-tf ) 
0 
0 
0 


<5 co  -  [6ti  ,  6 12  ,  <5 t^p 


] 


T 


(4.15) 


(4.16) 
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From  equation  (4.5),  (4.6)  and  (4.7) 


M>(tf)  = 


M. 
3X 


=  [0] 


Afi(tf)  9X 


tf  =  [0] 


3^1 

=  3X“ 


t  =  [0  0  0  -1  0] 
Lf 


X^2  ( ) 


dip, 


3X 


.  =  [.06708  -.02913  .15437  0  0] 

T 


T  T  4>(tf)  j 

Xcj)fi(t)  =  A4>  ( t )  -  y  Xft(t) 


or  X4>^(t)  =  [0] 


T  T  Mtf)  j 

Xljjft(t)  =  Xl^(t)  -  Xft(t) 


T  T  «  T 

or  Xtf.fi ( t)  =  Xif(t)  +  Xfi( t)  (since  fl(t^)  =  -1) 


or  Xijifl(t)  =  xj(t) 


(4.17) 

(4.13) 

(4.19) 

(4.20) 

(4.21) 

(4.22) 

(4.23) 

(4.24) 

(4.25) 


Substituting  equation  (4.17)  through  (4.25)  into  equation  (A4-1.15), 
(A4-1.16)  and  (A4-1.17)  gives 


I'M  -[0] 
I  iff.  =[o] 


(4.26) 

(4.27) 
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3 

m  =  l 

S  =  1 


A* 


-1  .T  .  _di 


J 


s  T 

xipn(t)  L$  w  1 

:s-l 


(4.28) 


Equations  (4.26),  (4.27)  and  (4.28)  substituted  in  equation 
(A4-l.ll)  gives  the  final  variational  equation  for  the  control  as 

3  -1  T  -1 

<5oj(t)  =  J  W  L  (t)  Xif)fi(t)  dijj  (4.29) 

s=l  s 

where  dip  =  [dij^  j  d^]  (4.30) 

Equation  (4.29)  is  considerably  simpler  than  the  original  expression 
given  in  equation  (A4-l.ll).  Expressions  similar  to  equation  (4.29)  were 

f4l 

obtained  by  other  workers  alsoL  . 

ii)  The  nonlinear  model 

The  procedure  for  the  nonlinear  model  is  similar  to  that  for  the 

linear  one.  Here  F^"  is  given  as 

x  3 


11 

a21 

a31 (1+x4) 

a41X3 

0 

12 

a22 

^22  0 +x4 ) 

a42x3 

0 

i2(1+x4) 

^23^  "*”x4 ) 

a33 

( a4] x] +a42x2^ 

0 

1 3X3 

a23x3 

(a31 xl+a32x2 

0 

wo 

,  r-COSXr- 

1 5  5 

a25COSX5 

"a35Sinx5 

0 

0 

(4.31) 
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where  the  parameters  a^j  are  as  given  in  section  A2-1  of  appendix  A2 


L  = 


1  Od-j a( t-t-j )  - 1  Od i  A ( t- 1^ )  1  Od-j  A ( t- 1^ )  -5d-|  A(t-t.p) 

1 Od^A ( t- ti )  - 1 Od^A ( t- t^ )  1 Od^A ( t- t^ )  - 5d^A ( t- t^) 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


(4.32) 


where  the  nominal  control  sequence  is  taken  as 


uQ(t)  =  5 
uQ( t)  =  -5 
uQ( t)  =  5 
uQ ( t)  =  -5 


t-|  <  t  <  t2 

<  t  <  tg 

tg  <  t  <  tf 


The  terminal  constraints  are 


=  "x4  =  0 

^  =  .0936667  x1  x3  -  .0406667  x2  x3  =  0 

Xij^  =  [0  0  0  -1  Of 


(4.33) 


(4.34) 

(4.35) 

(4.36) 


ANGLE  (deg  ) 
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X*2  =  [.0936667  x3(tf),  -.0406667  x3(tf),  (.0936667  x.,(tf)  - 

.0406667  x2(tf)),  0,  0]T  (4.37) 

The  rest  of  the  procedures  follow  that  for  the  linear  system 

4.3  Results  comparing  steepest  descent  with  quasi -optimal  schemes 

For  a  10%  torque  step  (figures  4.2  and  4.3)  it  can  be  seen  that 
there  is  virtually  no  difference  in  the  response  between  the  closed  loop 
scheme  and  the  optimal  scheme  (by  S.D.  method)  except  near  the  point  where 
the  control  switches. 


Fig.  4.2  Angle  time  characteristics  for  a  10%  torque  step 

(linearized  system). 


< 


t 


ACCLN.  (xio 
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Fig.  4.3  Phase  plane  corresponding  to  Fig.  4.2. 

1.  Optimal  control. 

2.  Closed  loop  control . 

There  is  no  doubt  that  the  closed  loop  (C.L.)  scheme  gives  optimal 
control  when  L(t)  is  zero.  For  small  disturbances,  the  magnitude  of 
the  term  L(t)  is  sufficiently  small  compared  to  the  input  term  to  yield 
an  almost  optimal  control  for  a  piecewise  constant  approximation. 
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Table  4.1  shows  that  as  the  integration  step  size  is  reduced,  the  time 
required  for  stabilization  is  less.  It  is  obvious  that  the  closed  loop 
scheme  depends  on  the  integration  step  size  and  does  not  foresee  the 
switch  curve  so  that  it  is  quite  likely  that  the  control  switches  several 
times  near  the  correct  switch  point  (figure  4.3).  This  can  be  avoided 
by  taking  larger  step  sizes,  in  which  case  the  'switch  miss'  may  be  quite 
large  resulting  in  at  least  another  extra  switch.  The  results  of  Table 
4.1  and  4.2  show  this.  The  best  way  to  avoid  the  intermediate  switches  is 
to  provide  a  'deadzone'  so  that  the  control  does  not  switch  if  the 
magnitude  of  the  switch  function  l  becomes  smaller  than  a  certain  pre¬ 
determined  quantity  e. 


Scheme 

Step  size 

(sec) 

No.  of 

swi tches 

Total  cost 

(sec) 

%  Error 

Remarks 

Steepest 
descent  (S.D.) 

.0005 

1 

.3415 

- 

- 

Closed-loop 

(C.L.) 

.0001 

1 

.3427 

.352 

Intermedi ate 
switches  (I.S) 

Closed-loop 

.0002 

1 

.3428 

.382 

Fewer  I.S. 

Closed-loop 

.0005 

1 

.343 

.44 

Still  fewer 

I.S. 

Cl osed-loop 

.001 

2 

.366 

6.7 

Almost  no 

I.S. 

Table  4.1  10 l  torque  step,  linearized  system  (ceiling  voltage 

±3  p.u . ) 


ANGLE  (deg) 
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For  large  disturbances,  the  results  do  not  agree  very  closely  to 
the  optimal  ones  as  shown  in  figures  4.4,  4.5  and  4.6,  which  are  for 
a  30%  input  torque  step.  This  is  approximately  the  maximum  load  the 
machine  can  deliver. 


Fig.  4.4.  Angle  time  characteristics,  30%  torque  step 

(nonl inear  model ) . 

1.  Optimal  control. 

2.  Cl osed  loop  control . 


* 
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Fig.  4.5  Speed  deviation  vs  time  corresponding  to  Fig.  4.4. 

1.  Optimal  control. 

2.  Closed  loop  control . 


During  and  immediately  after  a  large  disturbance  on  the  system, 
the  constraint,  that | L(t)/b( t) | <1 ,  may  not  be  met.  It  was  observed  in 
reference  [5]  that  if  for  a  small  fraction  of  the  total  time  under 
consideration,  L(t)  does  not  satisfy  the  constraint,  the  scheme  will 
yield  a  sub-optimal  solution.  However,  if  the  constraint  is  not  met  for 


(i-0 1  x )  ‘N1DDV 
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0.60  -r 


Fig.  4.6  Phase  plane  plot  corresponding  to  Fig.  4.4. 

1.  Optimal  control. 

2.  Closed  loop  control . 


0.70 
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all  t€:[tQ,  t^],  the  scheme  fails  to  give  any  information  about  the 
optimal  control.  Another  factor  to  be  considered  here  is  that  L(t) 
cannot  be  considered  as  a  disturbance  term  as  in  reference  [5]  since  it 
contains  output  feedback  s  .  In  such  a  case  the  solution  will  be  close 
to  optimal  only  if  the  dependence  of  L  on  6  is  small.  However,  the 
nonlinear  system  can  always  be  linearized  and  the  optimal  control  can  be 
applied  on  the  basis  of  a  linearized  system. 


Di sturbance 

Scheme 

Step  size 

No.  of 
swi tches 

Total  cost 
(sec) 

Error 

Remarks 

20%  torque  step 

S.D. 

.001 

1 

.33144 

- 

- 

20%  torque  step 

C.L. 

.001 

2 

.386 

16.3% 

* 

I.S. 

30%  torque  step 

S.D. 

.001 

1 

.409912 

- 

- 

30%  torque  step 

C.L. 

.001 

2 

.4788 

16.8% 

I.S. 

*  Refer  to  Table  4.1. 

Table  4.2  Nonlinear  system  (ceiling  voltage  ±5  p.u.). 


4.4  Effect  of  time  delay 

So  far  it  has  been  considered  that  the  control  can  be  instantly 
applied  to  the  system  immediately  after  it  is  found  either  by  optimal 
or  quasi-optimal  schemes.  In  general,  the  switching  times  of  the  control 
are  changed  from  the  nominal  ones,  by  the  delay  introduced  by  computation 
(includes  instrumentation  delay),  by  the  actuator  (includes  exciter)  and 
by  the  variation  of  the  parameters  of  the  system.  The  first  two  delay 


. 

- 

• 


* 
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the  switching  of  the  control  and  the  last  one  might  cause  the  switching 
to  occur  ahead  of  the  nominal  value  or  after  it.  However,  it  is 
considered  that  the  effect  of  variation  of  parameters  on  the  control  is 
small  compared  to  the  others.  In  this  section  the  effect  of  time  delays 
on  the  stability  limits  of  a  synchronous  generator  is  investigated. 


r 


~i 


Fig.  4.7  The  system  block  diagram. 

Two  different  cases  are  considered  here  -  the  effect  of  time  delays 
with  fixed  excitation  strategy  (or  the  so  called  open  loop  optimal  control); 
and  the  effect  on  closed  loop  excitation  control. 
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i)  Fixed  strategy 

The  switch  curve  for  a  particular  disturbance  is  found  by  the 
closed  loop  scheme.  The  switch  curve  crosses  the  time  axis  at  several 
points  and  these  are  stored  (this  is  equivalent  to  storing  the  control). 

It  is  also  necessary  to  store  the  switch  points  in  the  'deadzone'  (where 
the  magnitude  of  the  switching  function  is  less  than  10"^)  since 
observation  for  a  longer  period  of  time  is  needed  when  the  control  is 
delayed. 

The  second  step  is  to  delay  the  control  by  x  seconds  from  the 
predetermined  control  strategy  making  sure  that  the  field  voltage  remains 
normal  for  t  <  xsec.  Different  values  of  x  have  been  considered  for  the 
nonlinear  machine  model  with  20%  torque  step.  It  is  observed  that  x  =  .215 
seconds  is  in  the  neighbourhood  of  the  critical  delay  when  the  system 
appears  to  be  in  a  limit  cycle.  For  x  greater  than  .25  seconds,  the 
system  is  unstable.  The  phase  plane  plots  are  shown  in  figures  4.8,  4.9, 
4.10  and  4.11  for  different  values  of  x. 

The  system  without  any  time  lag  returns  to  the  equilibrium  point 
in  about  .4  seconds.  The  field  voltage  after  that  simply  bangs  around 
in  the  'deadzone'  with  no  net  effect.  The  disturbance  considered  is 
small  and  in  the  absence  of  any  stabilizing  signal,  the  machine  may  remain 
stable  up  to  .92  seconds.  Due  to  these  reasons  a  relatively  large  time 
does  not  seriously  deteriorate  the  stability  of  the  system  for  the  fixed 
excitation  strategy  which  has  little  importance  for  on-line  applications. 


. 

' 


Fig.  4.8  Velocity  vs  acceleration  plot  for  no  time  lag  (t  =  0) 


ACCLN.  (  p.u. ) 
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0.04  x 


Fig.  4.10  Phase  plane  plot  for  t  =  .215  Sec. 
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Fig.  4.11  Phase  plane  plot  for  x  =  .25  Sec. 
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ii)  Closed  loop  strategy 

The  quasi -optimal  control  based  on  the  system  states  is  computed 
and  stored  at  each  point  of  integration.  The  control  is  applied  to  the 
system  after  x  seconds.  The  delayed  control  results  in  new  system  states 
and  the  control  in  the  subsequent  steps  is  found  on  that  basis.  Since 
the  control  depends  on  the  measurement  of  actual  states,  it  is  expected 
that  the  time  lag  allowable  is  less  in  this  case.  For  the  same 
disturbance  as  in  the  previous  case  (open  loop  case),  a  delay  of  about 
50  ms  is  critical.  The  system  becomes  unstable  for  a  time  lag  of  60  ms. 
For  modern  excitation  systems  and  measurement  equipment,  this  is  an 
appreciable  time  from  a  stability  viewpoint.  The  phase  plane  plots  are 
given  in  figures  4.12  and  4.13  for  different  values  of  time  lags. 


Fig.  4.12  Velocity  vs  acceleration  plot  for  x  = 


.05  Sec. 
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The  last  case  considered  is  that  no  control  was  applied  up  to  about 
.92  seconds  followed  by  a  no-lag  bang  bang  control.  The  system  required 
about  2.5  seconds  to  stabilize.  This  shows  that  a  longer  initial  delay 
is  allowable  even  with  the  closed  loop  scheme.  The  critical  initial  delay 
is  seen  to  be  about  .95  seconds. 


* 


* 


107 


1. 


2. 


3. 


4. 


5. 


6. 


7. 


8. 


REFERENCES 

A.E.  Bryson  and  W.F.  Denham,  "A  Steepest  Ascent  Method  for 
Solving  Optimum  Programming  Problems",  Trans.  ASME, 

Ser  E,  Vol.  29,  J.  Appl .  Mech.,  pp.  247-257,  1962. 

W.F.  Denham  and  A.E.  Bryson,  "Optimal  Programming  Problems 

with  Inequality  Constraints  II:  Solution  by  Steepest 
Descent",  AIAA  J.,  Vol.  2,  pp.  25-34,  1964. 

R. F.  Vachino,  "Steepest  Descent  with  Inequality  Constraints 

on  the  Control  Variable",  J.  SIAM  Control,  Vol.  4,  No.  1, 
pp.  245-261,  1966. 

S. M.  Miniesy  and  E.V.  Bohn,  "Optimum  Network  Switching  in 

Power  Systems",  IEEE  Trans.  Power  App.  and  Systems, 

Vol.  PAS-90,  No.  5,  pp.  2118-2123,  Sept/Oct  1971. 

R.  Oldenburger,  "Optimal  Control",  Holt,  Rinehart  and  Winston, 
1966. 

A.E.  Bryson  et  al.,  "Optimal  Programming  Problems  with  Inequality 
Constraints  I:  Necessary  Conditions  for  Extremal 
Solutions",  AIAA  J.,  Vol.  1,  pp.  2544-2550,  1963. 

R. M.  Van  Slyke,  "Mathematical  Programming  and  Optimal  Control", 

Ph.D.  Thesis,  University  of  California,  Berkeley,  1965. 

S.  Dreyfus,  "Variational  Problems  with  State  Variable  Inequality 

Constraints",  RAND  Corporation,  1956. 


' 


108 


APPENDIX  A4-1 

METHOD  OF  STEEPEST  DESCENT 

Many  computer  procedures  for  solving  the  problem  stated  in 
Section  (3.2)  are  available  in  the  1  i terature^1  ’2 ,3 ,6^ .  Van  Slyke^ 
in  his  Ph.D.  thesis  made  a  detailed  study  of  mathematical  programming 
techniques  for  optimal  control  problems.  These  are  not  presented  here. 

Both  Pontryagin's  minimum  principle  (P.M.P)  and  the  method  of 
calculus  of  variations  can  be  applied  to  solve  the  problem  stated  in  the 
previous  section.  These  are  indirect  methods  in  that  the  control  is  a 
function  of  some  multipliers, values  of  which  are  not  known  in  general. 
P.M.P  uses  a  set  of  adjoint  variables  initial  values  of  which  are 
generally  not  known.  Newton  Raphson  method  may  be  used  to  search  these 
initial  values.  Convergence  is  fast  but  the  guess  for  the  initial 
costate  variables  should  be  quite  close  to  the  optimal  ones.  This  method 
is  seriously  handicapped  in  that  it  is  usually  difficult  to  make  a  good 
guess  for  the  initial  adjoint  variables. 

Iterative  methods  such  as  steepest  descent  and  conjugate 
gradient  methods  based  on  the  calculus  of  variations  deal  directly  with 
the  control  and  yield  the  optimal  solution  in  the  limit.  One  disadvantage 
with  these  methods  is  that  a  local  maximum  (or  minimum)  solution  may  be 
reached  instead  of  the  global  one.  The  conjugate  gradient  method  is 
faster  in  convergence,  but  needs  an  additional  set  of  influence  functions 
if  the  final  values  of  the  states  are  not  free. 


. 


- 


* 
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Due  to  the  simplicity  of  the  steepest  descent  method,  the  problem 
stated  is  solved  with  this  method.  It  can  be  seen  that  the  right  hand 
side  of  equation  (3.1)  is  linear  in  the  control  variable  u(t).  The 
optimal  control  for  such  a  problem  is  bang  bang  by  P.M.P.  The  only 
unknown  for  these  problems  is  the  optimum  value  of  the  switching  times 
of  the  control  u(t).  The  steepest  descent  algorithms  constructed 
by  Bryson  and  Denham^  and  later  modified  by  Vachino^  to  take  into 
consideration  control  problems  with  inequality  constraints,  have  been 
used  here.  The  problem  is  formulated  as  follows: 


Choose  the  control  function  u(t)  from  a  class  of  piecewise 


continuous  functions  of  time  such  that  |u(t)|<  1  which  takes  the  system 
described  by  the  vector  valued  differential  equation 


X  -  F(X,  u)  =  0 


(A4-1.1) 


from  its  initial  state  X_(o)  at  time  t  ,  to  its  intended  final  state, 
such  that  it  satisfies  the  vector  valued  terminal  condition 


^(X(tf),  tf)  =  0 


(A4-1.2) 


and  minimizes  the  cost  index  (Mayer  formulation) 


4>(X(tf),  tf) 


(A4-1.3) 


The  time  t^  is  chosen  as  the  first  time  that  one  of  the  terminal 
conditions,  referred  to  as  a  stopping  condition, 


fi(^(tf ) ,  tf )  =  0 


(A4-1.4) 


is  satisfied. 


* 


no 


Using  the  Heaviside  step  function^^ ,  the  equation  of 
for  equation  (A4-1.1)  can  be  written  as 


where 


6x(t)  =  F  6x  +  L5a) 

A 


L  _[(Fs(i+l)  "  Fsi 


s  =  1,2, 


N 


i  =  1 ,2,  . n 

N  is  the  total  no  of  distinct  discontinuities  in  u(t) 
for  the  system  of  order  n 
6oj  =  [ 6 1 1  ,  6t2,  . 6tN^T 

a( t-t  )  is  a  delta  function  occuring  at  the  switching 
instants  t  =  t 


The  influence  functions  are  given  by 


_T 

X  =  -  Fx  X 


Choosing  the  terminal  conditions  on  the  adjoint  vector 


V1^ 


W 


9X 

djj_ 

9X 

9X 


the  variation  of  control  w  (the  switch  times)  in  the  final  form 


variati on 

(A4-1.5) 

(A4-1.6) 


(A4-1.7) 

(A4-1.8) 

(A4-1.9) 

(A4-1.10) 


obtained  as 


- 


. 


m 


N+l 

6co(t)  =  ±  J  W 
S  =  1 


-1 


LI(t)CV(t)  -  V(t)  v1  lJ 


;dp)2  V  i^-1  d* 

-  iT  i  -1  i 

4><f)  ip<p  ipijj  ip( j> 


where 


j 


N+l 


+  l  W'1  lJ(t)  X  (t)  I  ~ 1 


-1 


s=l 


ipfi  tyip 


T,  x  T  i(tf)  T 


T  T  ^  (  "tf  )  T 

V'1  ■*»">-  SP^"’ 


rt, 


(dp) 


6  0)(  T  )  W  6  w(  T  )  dx 


J 

N+l 

^  '  S  =  1 


r  s 


J 


xl„  L  W'1  lT  xAn  dT 


s  - 1 


t 

N+l  f  s 

‘  J, 


XTn  L  W'1  Li  XAn  dx 

ipffi  S  S  (f>fi 


J 


's-1 


N+l  (  s 

I 

s=l 


x|0  L  W'1  lT  Xlo  dx 


•s-1 


(A4-l.ll) 

(A4-1.12) 

(A4-1.13) 

(A4-1 .14) 

(A4-1.15) 

(A4-1.16) 

(A4-1.17) 


Computing  Procedure 

a)  Compute  the  nominal  path  by  integrating  the  differential 
equations  with  a  nominal  control  variable  program  and  fixed  initial 
conditions.  Store  the  solution. 

b)  Compute  the  X^  functions  all  at  the  same  time  by  integrating 
the  adjoint  differential  equations  backwards,  evaluating  the  partial 


- 
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derivatives  on  the  nominal  path  by  reference  to  (a). 

c)  Simultaneously  with  (b)  calculate  the  quantity  x ^  and 
perform  integrations  (backwards)  leading  to  the  number  I  . 

d)  Select  desired  terminal  conditions  d^  to  bring  the  next 
solution  closer  to  \p  =  0  than  were  achieved  by  nominal  path. 

e)  Using  the  values  of  d^  chosen,  so  that  the  final  time  does 
not  exceed  the  predetermined  value,  calculate  5(u(t).  If  t^  exceeds  the 
predetermined  value  go  over  again  to  (d). 

f)  Obtain  a  new  nominal  path  by  using  w  =  u  ^  and 
repeat  process  (a)  through  (f)  until  terminal  conditions  \j>  =  0  are 
satisfied . 

j 


CHAPTER  5 


AM  EXTENDED  CASE 


This  chapter  considers  an  extension  of  the  method  developed  in 
chapter  3  to  a  two  machine  system.  The  closed  loop  control  scheme  for 
the  two  machine  system  is  developed  in  section  1.  The  results  obtained 
by  such  control,  for  a  number  of  cases,  is  discussed  in  section  2. 

5.1  Closed  loop  control  for  two  machine  system 

For  stable  operation  of  synchronous  generators  in  a  power  system, 
the  relative  velocity  and  acceleration  between  any  two  machines  in  the 
system  should  drop  to  zero  following  a  disturbance  while  the  relative 
rotor  angular  positions  should  remain  between  0  and  tt/2  radians.  So 
for  two  machine  system,  the  optimal  control  problem  for  stability  can 
be  formulated  as  follows: 

Given  the  initial  states  Xjo)  ;  m  =  1 ,2  and  the  following 
information  on  the  target  states 

n-j  (t.p)  -  n 2 ( t^r )  -  0 

(tf )  -  ft2(tf)  =  0  (5.1) 

0  <  <5-| (tf )  -  «2(tf)  <  tt/2 

where  nm  =  x4m  and  6m  =  “o  x5m 


. 
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Find  the  admissible  control  u  (t)  which  forces  system  (2.22)  from  the 
given  initial  states  to  the  desired  final  states  in  the  minimum  possible 
time. 

The  procedure  for  the  two  machine  system  is  similar  to  that  for 
the  single  machine  case  and  is  briefly  outlined  as  follows: 

For  each  machine,  differentiate  the  swing  equations  (2.23)  with 
respect  to  time  and  substitute  the  voltage  current  relations.  Then  the 
following  set  of  differential  equations  is  obtained 


p3(yv  =  M1) +  bi(t)  yo 


(5.2) 


P3(62/“0)  =  L2(t)  +  b2(t)  u2(t) 


(5.3) 


Subtracting  equation  (5.3)  from  (5.2)  yields 


P3(6i/.0)  -  p3(62/co0)  =  L-,(t)  -  L2(t)  +  b-j  (t)  u-j(t) 

-b2(t)  u2(t) 


(5.4) 


or 


p^  ( 6/ ajQ )  =  L(t)  +  u(t) 


(5.5) 


where  6  =  6 -j  -  6^ 


L ( t )  =  L-j  (t)  -  L2(t) 
u ( t )  =  b-j  ( t)  u -j  ( t )  -  b2(t)  u 2 ( t ) 


(5.6) 


As  before,  assuming  x-j  =  5/w0,  equation  (5.5)  gives  the  set  of  differential 
equations 
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X1  =  x2 

x2  =  x3 
x3  =  u0(t) 


where  uQ(t)  =  L1(o)  -  L2(o)  +  b1(o)  u,  (t)  -  b2(o)  u2(t) 
and  uomin  i  “0(t)  -  uomax 


(5.7) 


(5.8) 


Defining  the  following  quantities  as 

X-2 

1  =  x? 


2  ”  2u  •  *  lf  X3>0 

omm 

x32 

l  "  *2  -  2\j  ’  if  X3<0 


h  -  X1  ' 


omax 


2n  x2x3  ,  1  ^  ;(  3 


a 


3a 


2  3 


(5.9) 


(5.10) 


(5.11) 


X2X3  .1.3 


^  =  X1  '  “  +  3a2  *3 


(5.12) 


where  a  =  L(o)  +  uQmin  if  l  >  0 

a  =  L(o)  +  u  if  V  <  0 

a  v  ;  omax  L 


(5.13) 

(5.14) 


Then  the  control  scheme  is 

uo(t)  ”  uomax 

if  ^  <  0 

(5.15) 

uo^)  "  uomin 

if  l2  >  0 

(5.16) 

Otherwise,  uQ(t)  =  uQmax 

if  I  <  0 

(5.17) 

116 


(5.18) 


As  in  the  case  of  a  single  machine  system, the  constraint  on  L(t)  that 


L  ( t )  |  <  I  u  ( t ) 


(5.19) 


should  be  satisfied.  This  implies  that  u  .  <  0. 

orm  n 

Substituting  equation  (5.15)  through  (5.18)  the  initial  value  of  the 
quasi -optimal  control  u-j(t)  and  u ^ ( t )  is  obtained  from  equation  (5.8). 
The  process  is  continued  until  the  desired  final  states  are  reached. 

The  proportional  control  for  the  two  machine  case  is  obtained  in  a 
similar  way  as  for  the  single  machine  case. 

5.2  Discussion  of  results 

The  system  configuration  as  given  in  figure  2.3  is  considered. 

The  machines  are  interconnected  by  a  double  circuit  transmission  line 
and  feed  an  impedance  load  from  the  intermediate  bus.  The  two  machines 
considered  are  identical  but  carry  different  loads.  Machine  II  has  less 
load  than  machine  I.  Faults  and  disturbances  are  considered  on  machine  I 
only.  The  following  cases  are  considered.  (Data  for  the  two  machine 
system  are  given  in  section  A5-1  of  appendix  A5). 

i)  100%  torque  pulse  for  3  cycles. 

The  angle  time  characteristics  in  figure  5.1  shows  that  the 
application  of  bang  bang  control  results  in  a  near  deadbeat  response. 

This  changes  the  original  load  sharing  between  the  two  machines.  The 


. 
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proportional  control  based  on  the  relative  velocity  and  acceleration 
of  the  system  gives  a  heavily  damped  response,  bringing  the  system  to 
the  original  operating  point.  The  intermediate  bus  voltage  returns  to 
normal.  Though  the  relative  velocity  of  the  machines  becomes  zero,  the 
system  speed  up  by  about  .4%  (figure  5.4).  On  a  real  system,  this  will 
not  happen  since  the  governing  system  will  act  to  bring  down  the  speeds 
of  the  individual  machine.  Figures  5.2  and  5.3  give  the  phase  plane  plot 
and  figure  5.5  shows  the  variation  of  the  terminal  voltage  of  the 
individual  machines. 


1 


TIME  (sec) 

Fig.  5.1  Relative  angle  time  characteristics  for  100%  torque  pulse 

1.  Bang  bang  control. 

2.  Proportional  control. 


’ 


REL.  VELOCITY  (p.u.) 
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Fig.  5.2  Relative  velocity  vs.  relative  acceleration  plot 

corresponding  to  Fig.  5.1. 

1 .  Bang  bang  control . 


2.  Proportional  control. 


Fig.  5.3  Relative  angle  vs  relative  velocity  plot 

corresponding  to  Fig.  5.1. 

1 .  Bang  bang  control . 

2.  Proportional  control. 
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Fig.  5.4  Velocity  vs  time  with  proportional  control, 

100%  torque  pulse. 

1.  Machine  I 


2.  Machine  II 


TERMINAL  VOLTAGE 
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TIME  (sec) 


Fig.  5.5  Terminal  voltage  variation  with  proportional  control, 

100%  torque  pulse. 

1 .  Machi ne  I 


2.  Machine  II 


FIELD  VOLTAGE 
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Fig.  5.6  Field  voltage  characteristics  with  proportional  control, 

100%  torque  pulse. 

1 .  Machine  I 


2.  Machine  II 


- 
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ii)  150%  torque  pulse  for  3  cycles 

This  is  a  very  severe  disturbance.  The  bang  bang  control  results 
in  an  overshoot  and  an  undershoot.  The  intermediate  bus  (load  bus) 
voltage  is  very  much  reduced  resulting  in  a  totally  different  load  sharing 
between  the  machines.  The  proportional  control  results  in  a  heavily 
oscillatory  response  (fig.  5.7)  and  a  return  to  nearly  normal  load  sharing. 
Figure  5.8  gives  the  phase  plane  plot. 


Fig.  5.7  Relative  angle  time  characteristics  for  150%  torque  pulse. 

1 .  Bang  bang  control . 

2.  Proportional  control. 


REL.  VELOCITY  (  p.u.) 
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Fig.  5.8  Phase  plane  plot  corresponding  to  Fig.  5.7 

1 .  Bang  bang  control . 

2.  Proportional  control. 
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Hi)  Three  phase  short  circuit  at  the  transformer  end  of  the  first 
machine  cleared  in  3  cycles  (underexcited  case) 

The  second  machine  is  normally  excited.  Both  the  machines  are 
relatively  lightly  loaded.  (Machine  I  is  loaded  at  .73  per  unit). 

Figure  5.9  shows  that  both  bang  bang  and  proportional  control  results 
in  similar  type  of  oscillations.  The  proportional  control  restores  the 
load  bus  voltage  and  machine  internal  voltage,  while  the  bang  bang 
control  reduces  them.  Figure  5.10  shows  the  phase  plane  plot  and  figure 
5.11  shows  the  variation  of  the  individual  rotor  angles  with  proportional 
control  only.  It  can  be  seen  that  the  rotor  angles  continue  to  increase 
with  time;  this  is  due  to  the  off-set  in  velocity  of  the  individual 
machines.  Again,  this  will  not  occur  if  governor  action  is  taken  into 
consideration. 


TIME  (sec) 

Fig.  5.9  Relative  angle-time  characteristics  for  3  -  <j>  fault 

on  under  excited  machine. 

1 .  Bang  bang  control . 

2.  Proportional  control. 


. 


REL.  VELOCITY  (p.u.) 
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Fig.  5.10  Relative  angle  vs  velocity  plot  corresponding 

to  Fig.  5.9. 

1 .  Bang  bang  control . 

2.  Proportional  control. 


ANGLE  (cleg) 
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Fig.  5.11  Angle  time  characterise cs  with  proportional  control. 

1 .  Machine  I 

2.  Machine  II 

iv)  Three  phase  short  circuit  on  first  machine  (normal  excitation). 

In  this  case  machine  I  is  loaded  about  .85  per  unit  and  machine  II 
about  .4  per  unit.  Figure  5.12  shows  that  bang  bang  control  deteriorates 


» 
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the  machine  and  bus  voltages  so  much  that  when  the  system  tends  to 
stabilize  after  about  .8  seconds,  the  maximum  power  delivered  by  the 
two  machines(at  these  voltages)is  much  less  than  that  demanded  by 
the  load,  resulting  in  instability.  Proportional  control,  however, 
retains  the  machine  as  well  as  bus  voltage,  so  that  the  system 
stabilizes  in  about  1.6  seconds. 


Fig.  5.12  Relative  angle  time  characteristics  for  a  3  -  <f>  fault 

on  the  normally  excited  machine. 

1.  Bang  bang  control. 

2.  Proportional  control. 

v)  The  effect  of  local  control. 

A  100%  torque  pulse  for  3  cycles  is  again  considered  on  the 
first  machine.  Figure  5.13  shows  the  relative  angle  vs  time  characteristics. 


* 
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Curve  1  is  with  bang  bang  control  on  both  the  machines,  whereas  curve  2 
is  with  control  on  the  first  machine  only,  second  machine  with  normal 
field  voltage.  Curve  3  shows  the  response  when  control  is  applied  on 
machine  II  only.  It  can  be  seen  that  applying  control  on  the  first 
machine  only  takes  less  time  as  compared  to  that  when  control  is  applied 
on  both  of  them.  This  demonstrates  that  local  disturbances  are  better 
taken  care  of  by  local  control.  In  this  case  the  local  control  also 
depends  on  the  global  variables.  Figure  5.14  shows  the  phase  plane  plot 
corresponding  to  curves  2  and  3  in  figure  5.13. 


Fig.  5.13  Relative  angle-time  characteristics  for  100%  torque  pulse. 

1.  Control  on  both  the  machines. 

2.  Control  on  Machine  I  only. 

3.  Control  on  Machine  II  only. 
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Fig.  5.14  Phase  plane  plot  corresponding  to  Fig.  5.13. 

1.  Control  on  Machine  I  only. 

2.  Control  on  Machine  II  only. 
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APPENDIX  A5-1 

TWO  MACHINE  SYSTEM  DATA 


The  two  machines  considered  are  identical  and  their  parameters 
are  given  in  section  A3-3  of  appendix  A3.  Two  system  conditions  were 
considered  -  (1)  both  machines  normally  excited,  and  (2)  machine  I 
under  excited. 

The  operating  quantities  for  the  under  excited  case  is  given  as 


M/C  I  M/C  II 


nfd 

1.42348 

Vd 

1.77935 

jd 

.231041 

1d 

.0066098 

\ 

.987381 

\ 

.19621 

6 

60° 

6 

10° 

6t 

.7696 

6t 

.9983 

Po 

.72976 

Po 

.1957 

Q0 

-.27684 

Q0 

-.01021 

v  =  1 

Z.  =  .852  /_-.712  rad. 


Operating  quantities  for  the  normally  excited  case  are 


■ 
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M/C  I  M/C  II 


]fd 

2.3131 

Vd 

2.224 

jd 

.4464 

]'d 

.2407 

]q 

.745 

1#q 

.3027 

6 

o 

o 

■=3- 

6 

15° 

et 

1.044 

et 

1.0927 

P 

0 

.8841 

P 

0 

.3599 

Qo 

.2025 

Q0 

.2215 

v  =  1 

=  .8246  j_. 0142  rad. 


CHAPTER  6 


CONCLUSIONS  AND  SUGGESTIONS  FOR  FURTHER  RESEARCH 

6.1  Summary  and  Conclusions 

In  this  thesis,  a  quasi -optimal  excitation  control  for  power  system 
stabilization  is  found  directly  as  a  function  of  state  and  parameters  of 
the  system.  It  is  suited  to  online  controls  requiring  relatively  little 
hardware,  either  digital  or  analog,  as  it  is  not  necessary  to  solve  the 
system  differential  equations.  Based  on  the  closed  loop  scheme,  two  types 
of  excitation  control  have  been  investigated: 

'i)  Bang  bang  control  scheme  -  This  provides  dead  beat  response  for 
small  disturbances  (figures  3.5  and  3.6).  However,  for  large  disturbances 
there  may  be  an  uncontrollable  period  resulting  in  a  sub-optimal  solution. 

The  exciter  ceiling  voltage,  which  limits  the  control  signal,  has  a  direct 
bearing  on  both  the  uncontrollable  period  and  the  speed  of  response  during 
the  controllable  period.  Increasing  the  exciter  ceiling  provides  better 
damping  of  the  transient  (figure  3.5).  If  a  bang  bang  control  is  used, 
some  deadband  must  be  provided  to  allow  normal  voltage  regulator  action  to 
take  place.  This  also  can  take  care  of  any  switching  chatter. 

ii)  Proportional  control  -  In  the  case  of  a  large  disturbance,  the 
exciter  is  rapidly  driven  to  its  ceiling  and  the  control  action  is  almost 
bang  bang.  As  the  system  settles  down,  the  exciter  voltage  variations 
decrease  and  the  effects  of  voltage  regulator  action  become  dominant. 

Dynamic  stability  is  retained.  Although  the  transient  is  longer  with 
greater  excursions  in  rotor  angle,  the  transition  to  normal  voltage  regulator 
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action  is  smooth  and  does  not  require  additional  control  equipment.  The 
demand  on  the  excitation  system  is  lower.  The  regulator  time  constants 
may  be  taken  into  account.  The  results  obtained  with  these  two  types  of 
control  are  presented  in  Chapters  3  and  5. 

The  results  obtained  with  the  closed  loop  quasi -optimal  scheme  are 
compared  with  those  obtained  by  a  steepest  descent  method.  For  small 
disturbances,  the  closed  loop  scheme  is  found  to  yield  an  optimal  solution. 

It  may  be  noted  that  although  the  steepest  descent  technique  was  used  to 
evaluate  the  closed  loop  scheme,  the  amount  of  computation  required  for 
steepest  descent  methods  precludes  online  control.  This  may  give  worse 
results  than  the  closed  loop  scheme  for  other  than  the  design  disturbances. 
For  the  application  of  steepest  descent  methods,  the  closed  loop  scheme 
offers  computational  advantages  in  that  the  results  provide  a  good  set  of 
initial  switching  times  for  the  steepest  descent  algorithms,  so  that 
convergence  is  rapid.  The  comparison  of  closed  loop  scheme  with  the  steepest 
descent  method  is  given  in  Section  4.3  of  Chapter  4. 

The  effect  of  instrumentation  and  related  time  delays  is  not  as  severe 
for  predetermined  control  strategies  as  for  the  closed  loop  scheme.  This  is 
shown  in  Section  4.4  of  Chapter  4.  Since  the  sensing  and  measuring 
instruments  as  well  as  the  static  excitation  systems  have  very  small  time 
constants,  their  effect  is  less  appreciable.  Even  a  relatively  long  initial 
time  delay  followed  by  a  bang  bang  excitation  strategy  retains  system 
stabi li ty. 

For  multimachine  systems,  the  control  is  obtained  based  on  the 
information  from  all  machines  in  the  system.  For  local  disturbances,  local 
control  may  be  more  effective,  in  which  case  the  amount  of  computation  is 
also  minimized. 


6,2  Suggestions  for  further  research 


Further  areas  of  investigation  include  the  extension  of  these 
methods  to  include  governor  and  prime  mover  dynamics  and  consideration 
of  damper  windings.  Tests  on  real  machines  should  be  performed  to  compare 
the  results  with  computer  simulation.  Methods  for  determining  L ( t )  and  b(t) 
or  reasonable  approximation  of  these  quantities  with  a  minimum  of  comput¬ 
ation  and  equipment  should  be  investigated.  The  sensitivity  of  these 
quantities  to  parameter  variations  are  also  of  concern. 

Modern  power  systems  are  very  complex  and  the  possibility  of 
monitoring  the  optimal  control  from  a  central  computer  should  be  explored. 
Study  of  local  control  involving  the  local  variables  only,  needs  concern. 

Further  investigation  of  local  versus  central  control  should  be 
made  as,  if  local  control  is  satisfactory  from  a  stability  viewpoint,  it 
would  be  superior  in  terms  of  cost,  equipment  and  complexity  and  speed  of 
control . 
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